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ABSTRACT 


A  preliminary  design  of  a  coal  burning  25 , CCO 
H?  naval  ship  propulsion  plant  was  completed.  The  system  con¬ 
sists  of  a  combined  Closed  Brayton  cycle  and  bottoming  R aricine 
cycle.  Coal  combustion  capability  is  achieved  by  using  a 
Pressurised  Fluidised  Bed  (??2)  heater.  A  thermodynamic 
analysis  and  volume  estimate  was  conducted  for  various  Brayton 
cycle  compressor  pressure  ratios  (  7TC  ) .  : 

The  results  show  that  a  maximum  thermodynamic  1 

efficiency  (  )  is  achieved  in  the  range  of  Tc  =  5  to  6 .  j 

Total  volume  on  the  other  hand  has  a  minimum  around  =  3.5 
and  climbs  steeply  above  that.  A  trade-off  exists  between 
efficiency  and  volume  with  a  good  design  range  of  tc  =  -  to  1 

The  design  point  was  ultimately  chosen  at  7Tr  =  a. 5  with  a  H-r  ~ 

Ac.  '  ' 

The  design  of  a  modern  coal  burning  combined  cycle 
offering  high  efficiency  at  both  rated  and  part  load  is  within 
the  reach  of  present  technology.  Heat  exchangers  and  turbo- 
machinery  design  are  well  within  the  state  of  the  art  and  ?F3 
technology  is  expected  to  come  on  line  within  the  next  few  years. 
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A3S TRACT 

A  preliminary  design  of  a  coal  burning  25.CCO 
HP  naval  ship  propulsion  plant  was  completed.  The  system  con¬ 
sists  of  a  combined  Closed  Brayton  cycle  and  bottoming  RarJcir.e 
cycle.  Coal  combustion  capability  is  achieved  by  using  a 
Pressurised  Fluid ized  Bed  (PF3)  heater.  A  thermodynamic 
analysis  and  volume  estimate  was  conducted  for  various  Bray ten 
cycle  compressor  pressure  ratios  (  7TC  )  . 

The  results  show  that  a  maximum  thermodynamic 
efficiency  (  7*  )  is  achieved  in  the  range  of  Tc  =  5  6. 

Total  volume  on  the  other  hand  has  a  minimum  creund  770  =  3.5 
and  climbs  steeply  above  that.  A  trade-off  exists  between 
efficiency  and  volume  with  a  good  design  range  of  Uc  =  1  to 
The  design  point  was  ultimately  chosen  at  7TC  =  ^.5  with  a  ^  s 
.13 . 

The  design  of  a  modern  coal  burning  combined  cycle 
offering  high  efficiency  at  both  rated  and  part  load  is  within 
the  reach  of  present  technology.  Heat  exchangers  and  turbo- 
macninery  design  are  well  within  the  state  of  the  art  and  PFB 
technology  is  expected  to  come  on  line  within  the  next  few  year 
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INTRODUCTION 


1.  Purpose .  \  Increasing  prices  ar.d  decreasing  supplies  of 
fuel  oil  has  created  a  need  for  highly  efficient  ship  pro¬ 
pulsion  systems  capable  of  utilizing  fuels  other  than  oil. 

The  purpose  of  this  thesis  is  to  perform  a  preliminary  design 
of  a  combined  cycle  propulsion  plant  eo.uipped  with  a  multi¬ 
fuels  combustor  capable  of  burning  coal  or  oil.  <r~~ 

2.  System  Description.  The  system  investigated  is  a  closed 
regenerative  3raytcn  Cycle  combined  with  a  bottoming  RarJcine 
cycle,  corking  fluids  for  the  cycles  are  air  and  water 
respectively.  The  heater  for  the  Bray ton  cycle  is  a  pressur¬ 
ized  fluidized  bed  (??3).  The  heater  cycle  is  essentially  an 
open  regenerative  Brayton  cycle.  This  thesis  was  conducted 
considering  the  heater  only  in  its  coal  burning  mode  of  opera¬ 
tion.  Work  done  by  the  Braytcn  and  Rank ine  cycles  will  drive 

a  common  shaft  through  a  double  reduction  type  gear.  An  enthalpy- 
enthropy  diagram  is  shown  on  Figure  1 . 

3 .  Scoce . 

a.  Objectives.  The  objectives  of  this  thesis  are  to: 

(1)  Perform  a  thermodynamic  analysis  of  the  system 
for  various  closed  3rayton  cycle  pressure  ratios  utilizing 
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realistic  values  for  component  efficiencies  ar.d  losses. 

(2)  Estimate  major  component  volumes  and  thus  obtain 
a  total  volume  for  set  of  parameters  listed  in  (!'  above. 

(3)  Select  a  set  of  parameters  for  design  based 
upon  the  above  results. 

(-)  Discuss  the  technical  feasibility  and  implication 
of  the  design  as  a  means  for  naval  ship  propulsion. 


t .  tester.  writer  la . 

(1)  The  system  parameters  selected  for  design  were 
those  which  gave  the  lowest  total  volume  within  an  acceptable 
range  of  thermal  efficiency. 

(2)  An  acceptable  range  of  thermal  efficiency  is 
defined  as  t  1  f»  of  the  maximum  efficiency. 


o .  limitations  to  S  core  . 

(1)  ;»o  component  or 'system  testing  was  carried  out. 

(2)  The  system  was  designed  with  the  aid  of  a 
computer  utilising  performance  parameters  of  currently  avail¬ 
able  equipment. 

(3)  The  concept  of  pressurised  fluidised  bed  (?F3) 
combustion  has  not  been  proven  in  a  shipboard  environment. 

The  thesis  was  conducted  based  on  the  assumption  that  a  ?F3 
system  would  be  put  to  sea  in  the  near  future. 


a 
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THERMODYNAMIC  ANALYSIS 


1 .  The  Braytcn  Cycle 

a.  ’■•/or  king  cycle.  This  is  a  closed  Bray  ton  cycle 
utilising  air  as  the  working  fluid.  Air  was  chosen  since 
this  is  proposed  for  surface  ship  application  and  would  thus 
eliminate  the  need  for  carrying  a  reserve  supply  of  gas  if 
something  like  Helium  or  GO,  were  used. 

Advantages  of  a  closed  cycle  over  an  open  cycle  are 


the  following: 

(1)  Compressor  suction  pressure  is  no  longer  iim:  c 
to  ambient.  By  increasing  the  working  pressure  of  the  eye 

a  significant  reduction  in  engine  site  can  he  achieved  foi 
specified  output. 

(2)  Higher  working  pressures  also  increases  the  heat 
transfer  capacity  of  the  working  fluid  thus  resulting  in 
smaller,  more  effective  heat  exchangers. 

(3)  The  use  of  a  high  pressure  fluid  gives  the  added 
advantage  of  regulating  the  cycle  pressure  level  hy  varying 
the  inventory  of  gas  in  the  loop.  If  this  is  done  in  a  manne 
which  kept  the  maximum  temperature  and  speed  constant,  the 
turbomachinery  vector  diagrams  would  remain  essentially  ’un¬ 
changed,  The  result  is  little  or  no  change  in  compressor  or 
turbine  efficiencies  with  varying  load,  thus  maintaining  high 
part  load  efficiency. 

A  Temperature  -  entropy  (T-s)  diagram  of  the  working 
cycle  is  shown  in  Figure  2. 


o 


~vj  U\Ul 
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Figure  2. 
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2 
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compressor  suction,  cooler  dischar 
compressor  discharge,  regen  inlet 
regen  outlet  (cold),  heater  inlet 
heater  outlet,  turbine  inlet 
turbine  outlet,  regen  inlet  (hot) 
regen  outlet  (hot),  WH3  inlet 
.vH3  outlet,  cooler  inlet 


heat  input  from  heater  cycle 
heat  released  to  Rardcine  cycle 
heat  released  to  sea  water  cooler 
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The  analysis  was  conducted  by  assuming  values  for 


T,  .  T,  ,  7c  ,'su  comor  ess  or  oolytrooic  efficiency) 

I  *•’  9 

polytropic  efficiency),  Z,  (regenerator  effectiveness),  and 

*  V 

(total  pressure  drop). 

The  specific  heat  was  assumed  to  vary  with  temperature 
according  to  the  following  relation: 


cine 


:p(T)  =  .2i?5  -  (3.759x13*^) 


2  ?  _  •  ' 

_  1  C  *  ~  C.  ^  ^  «  ...  •*  -  - 

J  ...  -\-.iJ-X-u  / 


:  or 


T  in  °H 


O- 


3?  in  BT'J/lbn-  ?. 


?he 


»  * 


\?--T 

•"  -a 
• 


a 

o 


"oil-owing  are  defined: 

compressor  pressure  ratio  \?~/?.] 

turbine  pressure  ratio  (?-/r\) 

-  '  a»  total  pressure  drop  around  closed  Bray  ton  loop 

power  (net)  from  closed  3rayton  cycle 
mass  flow  of  working  fluid  in  Brayton  cycle 

temp  in  °R  at  point  i 

enthalpy  at  point  1 

gas  constant 

pressure  ats .  at  point  i 

specific  volume 
regenerator  effectiveness 


The  basic  procedure  was  to  determine  ™  ^  *5'  *6' 

•  • 

Q-n,  and  .v^  for  various  pressure  ratios  (  7T\)  . 
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r  > 


Determination ' of  T,  ; rsf .  Dixon): 


Small  stage  polytropic  efficiency  is 


n  p: 

,C'  ' 


mce  tor  an  isantrcpcc  process  -  u=c: 


uG3  ^ >*»  * u  - 


T" =  RT/p  gives 


•Vp)  =  5_  1  12 


hus , 


.<  = 


■ — - 


v  substitute 


•'  f>  '  „  3 


2ii?5  In:  -  2  3 *  ?59xl . 


-  1  >  V 


2.553x:o"c)( 


O  1 

-  •  ~>j  ^xi  v  ;  (  D  5 


s  p  /  

p - M>  ->c 

';<?) 


then  be  determined  iteratively. 


Determination  of 


In  a  manner  similar  to  the  compressor  above 


C  ^2_i-  ■  R  **>  A  (A) 


Determination  of  ?  and  r ^  : 

_ d. _ 2 

T-  and  ?/  car.  de  found  from  an  assuned  regenera 
J  ^ 

effectiveness  (Z^)  using  the  following  relation: 
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V  i\ 


b 


Heater  Cycle.  This  is  an  open  Bray ton  cy 
combustion  air  as  the  working  fluid.  The  heater  i 
a  pressurised  fluidized  bed  utilizing  coal  as  fuel 
diagram  of  the  heater  cycle  is  shown  in  Figure  3 • 


le  also 

3  3  13 


US 


Major  components  of  the  cycle  consist  of  a  heater,  re¬ 
cuperator  and  turbine  driven  compressor.  The  turbine  is 
powered  by  exhaust  gases  leaving  the  recuperator. 

The  analysis  was  conducted  by  assuming  values  of  T,  , 
T?h.  jLC?)»  and  i  (equivalence  ratio}.  Both  the 

exhaust  products  and  combustion  air  were  assumed  to  be  perfec 
gases  with  specific  heats  varying  in  accordance  with  the  re¬ 
lation  given  in  i.(a). 


The  following  are  defined: 


i  1  Q  i* 

- 

heater  ccmcressor  pressure  ratio 

C  « i  * 

heater  turbine  pressure  ratio  (P,,,,/?^) 

- 

total  pressure  drop  around  heater  loop 

a 

a 

— 

mass  flow  of  air  in  heater  cycle 

a, 

a 

- 

mass  flow  of  fuel  in  heater  cycle 

( r  / k ,  £ 

- 

stoichiometric  fuel  -  air  ratio  of  fuel 

LKV 

- 

Lower  Heating  Value  of  fuel  used 

m 

*  ih 

- 

absolute  temp  at  point  i^ 

h,  u 

- 

enthalpy  at  point  i. 

efficiency  of  the  heater  cycle 

i 

Vh 

- 

erec 

- 

recuperator  effectiveness 

used 


14 


r* 


ambient 


figure 


point  lh: 
2h: 
3h: 
^h: 
5'n: 
6hs 
7h: 


compressor  ir.iet 

compressor  outlet,  recup  inlet  ( 

recup  outlet  (cold),  heater  inle 

max  air  temp  in  heater 

heater  outlet,  recup  inlet  (hot) 

recup  outlet  (hot),  turbine  inlet 

turbine  exhaust 


Qg*  heat  released  to  working  cycle 


heat  in  3  s,  LHV 
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ci  O 


Assume  the  heater  cycle  to  be  a  control  volume  :.~Z  . 

as  shown  in  Figure 


m^LHV 


v 


ma  nlh 


m.,)  h-. 
a  c  n 


ire  . 


leaving  the  C.7.  must  be  ecual  to  entering  the 


working  cycle  thus, 


Q3  =  m3  (hj,,  -  h-,)  *(na  -  af)(h 


,i..  -  h ) 

^h  5*t 


Performing  a  heat  balance  on  the  C.7.  gives  us 


3  3  *  ^a  +  h7h  =  •*  m  h 


a  Ih. 


Therefore , 


n  „  ** 


■3a 


“  |  (  I  +  (^'0s  ?)  n7*  “  $  L  1  v] 
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on  a  power  balance  between 


the  turbine  and  conpressor 


W  =  :U  ( no  net  work) 


^a  v*h2 


( -i 


)  = 


a  n  ~  / 

i  -  a; 


r'* 


n 


-  *7  Vs 


but 


I  '  >1 , 


^  n  ( 

—  /—  \  '■  '  *A'V 

Ti'n  =  '7h  »X"W 


*  il1  *  J  McV 


thus , 


”  -i  ~ i '  ' ' 


-  /s/  \  A-?  pv  -  r,V  ,'-T' 

=  I  in  /Ajs  0J  G-  j7k  i/th;  'J 

-  /  ,  \  /+~i  —  —  (r  /aF  •  “)  (7  \VT,'^‘  _ 

=  !:  ^f/,0s?]  -.--J  .iT<y 


or 


Ljr.  1  ik  tV'<k 


-  Cj.  '''> 
f.'_  \  V"-  /  ■/-- 


!  +1  |  + 

J  L 


(%),'?]  C,, 


'7  k 


/  \  7  ^  -  r ,  /AP\  -  4 

[  i  +  v'Vi  G  Vf  ;7k  L  '  *\  f3/'  i 


—  (  ,y  ■  U»  Z"'  ■>  : 

'  -  / 

"c 
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3y  substitution 


Thus  ^c, 
quantities 


found  through  known  or  assur.ec 


Te termination  of 


_  -r0K)  _ 


rec 


r 


2k  ~  tk 


where 


CN  =  t 


C,  = 


**K  a’ra*,n»*Q 
w.c  - —  i 


^  -  ‘-A 


+  [/  +(7*)$  ?\  ~S>K  'T^j 


C,.  and  T/.  can  be  found  once  'lr.  is  determined.  Assume 
ch  on  -p- 


that 


<i 


h  ‘3 


♦  1  00  Pi 
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Major  components  of  the  cycle  are  the  WK3,  condenser, 
steam  turbine,  and  feed  pump.  It  is  assumed  chat  me  feed 
pump  is  driven  with  power  from  the  turbine. 


The  analysis  was  conducted  by  assuming  values  of  T«  , 

1  T ,  { pump  iser.tropic  efficiency),  y-r'(i5;  (turbine  isen- 

tropic  efficiency),  and  ?  V.:H3  (waste  heat  boiler  efficiency 


The  following  are  defined : 


is 


pp 


or 


X 


p  _ 

is 


Y). 

IT 


cu 

o 


enthalpy  at  point  is 

absolute  temp  at  point  is 

pinch  point  on  3raytor.  side 
saturation  point  c n  steam  side 
mass  flow  ratio  (cu/m~) 

w  £ 

quality  of  steam  exiting  turbine 
absolute  pressure  at  point  is 

power  (net)  from  the  steam  cycle 

thermal  efficiency  of  the  combined  cycle 
mass  flow  of  steam  in  Rankine  cycle 


The 

(1) 

(2) 


Rankine  cycle  parameters  were  chosen  for  maximum  work, 
procedure  was  as  follows: 

Compute  T  ^  from  the  relation 


X3S  ‘  ‘6 

Assume  a  value  for  TOQ. 


-  AT 

T.  <  T  <  662.0°? 

1  pp 
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On 


(12)  Calculate 


from  the  relat 


•  *  *  I 


VCLUN'E  analysis 


1  .  Turbonachir.ary .  A  n  analytical  prcced  ure  for  the  minimum 
volume  design  of  turbomachinery  was  developed  by  lee  {ref. 
10) .  The  following  equations  and  mathematical  expressions 
summarize  that  portion  of  the  procedure  which  are  applicable 
to  compressor  and  turbine  designs  in  this  thesis.  Volume 
obtained  by  lee’s  method  neglects  the  volume  of  casing,  ir.le 
and  outlet  plenums,  and  bearing  assemclys. 

a.  Turbine.  Tree  vertex  blading  and  axial  exit  stages 
were  assumed  throughout  the  analysis. 


The  following  acre  defined: 


3t 

A 

>^■3 


'/<K 

rT 


n 


K 


a 


factor  of  safety 
max.  allowable  tip  speed 
flow  coefficient  at  tip 
hub-tip  ratio  at  inlet 
hub-tip  ratio  at  exit 
clearance  factor 
specific  volume  at  inlet 
max  tip  radius 
polytropic  coefficient 
iser.tropic  exponent 
density  of  blade  material 
stress  concentration  factor  at 
amplification  factor 
correction  factor  for  taper 


root 
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-p  =  r  r  L  - 


r*  m  =  P 

— >  i  *  u  b  i 


/  j  *  ;  r  J  ,'  '  7'  -•*'/,-  1 


3: 


-  -7 


0  / 


The  axial  chord  for  a  stags  is  given  by  3-  .  For  a 
large  number  of  stages,  summation  of  c-  for  the  who  Is  t; 
can  be  approximated  by  an  integral: 


'  V  ^ 


I  ^  V  ck 


\ 

f  ■ 
=< 


wnere 


-  K 


j n  \  .  _  /V_,  '  '  nN  o  OD  *N  i  /-pi  7 

^ l3'n y  "  1  \  •>7'  .  p'  j 


This  integral  can  be  evaluated  numerically  thus, 


\/„  ! 
yc  >r 


-  Z7T  r.'c 


1  It 


b.  Conoressor .  An  approach  similar  to  that  of  the 
turbine  was  used  for  the  axial  compressor.  Axial  exit  stages 
were  assumed  in  the  analysis. 

The  volume  was  computed  by  assuming  values  for  SF,  'J*., 

v 

'3-wiXo  CT  .  X,  J  .  *  ,  a,  x.  ,  Kt.  ,  9,  and  T„  iblade  material 

v.  u  (  m  -  y 


yield  stress) . 
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A 3  with  the  turbine  define 


;  > 


^  S°*V- 


where 


C\ 


.  i.  *'  '/■  \ 
“  z  \  '3J 


-•/. 


-  3JD  ‘  N, 


This  integral  can  be  evaluated  numerically  thus 


7  o  1 . 


•o  ^ 

U.  i  t  _ 


2.  Heat  exchangers .  The  volumes  of  the  Regenerator,  Re¬ 
cuperator,  Cooler,  and  Heater  were  computed  using  the  method 
of  minimum  volume  design  developed  by  lee  (ref.  1C).  In 
essence  lee  found  that  the  optimal  diameter  ratio  for 

these  heat  exchangers  can  be  established  as  a  simple  relation 

to  the  pressure  ratio  (  77  )  of  the  cycle.  The  results  of  the 

A  ° 

optimally  distributed  pressure  drops  (  ~=r  )  as  reported  by 
Lee  were  also  studied.  Parametric  equations  relating  pressure 
drop  to  pressure  ratio  were  developed  and  used  in  this  analyse 
No  consideration  has  been  made  on  computing  the  volume  of  the 
shell,  headers,  or  other  appendages.  Lee’s  relations  were 
developed  from  basic  heat  transfer  and  pressure  drop  relation¬ 
ships. 
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a.  Regenerator .  The  regenerator  is  assumed  to  be  a 
counter  flew  shell  and  tube  heat  exchanger.  High  pressure 
is  inside  the  tubes  as  is  the  case  for  current  design  crac 
tices . 


The  following  are  defined: 


_4?\ 

(  -t>  /  \;Q 


fi  pressure  drop  across  the  regenerator 
(both  steams) 

diameter  ratio  for  regenerator 

enthalphy  difference  across  the  interior 
stream  of  Reg. 

tube  inside  diameter 

log  mean  temp,  difference  of  hot  and  cold 
streams 

Prandlt  number 
density 

dynamic  viscosity 
thermal  conductivity 


A  P  \  'IX  v 

Values  for  "TyRo?  and  r£~y  Reg  can  be  estimated  from 
She  following  relations? 


=  i -5c  _  02q4TT  (for  — )  = 

(  pyReg  ,u  ic  Uor.^p/r 


.1,  ZR  =  .88) 


R-)r.S  -  -53?  *  .°9X  » e 
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Values  f:r:-r  and ' —  , _ can  be  estimates 

VP/P.ec  v-/^ec 

'allowing  relations 


/  a..3  \  _  ,'ap  \  f  a?\  /  ap  \ 

V  ■=»  A«".  ~  i  *  ~  v  7>  yr.Trj.  "  l,  ~7“ 


,£k_\ 

\  J> 


Subscripts  1  and  2  refer  to  quantities  inside  and  outside 
the  tubes  respectively .  A-  „  is  computed  as  follows: 

•  •*  «<  p  /**  * 


A  *  {%)  ~r  -r 


where  T  and  Y  are  as  previously  defined 


Then  compute  the  volume 


731  =  q  .**T>  -0>1 

‘‘■Rec  Sec  \  ?  yRec 


c.  ^oolar .  The  cooler  is  assumed  to  be  a  shell  and  tube 
heat  exchanger  with  the  gas  inside  the  tubes  and  salt  water 
outside . 


The  following  are  defined: 


(AL\ 

C  pycir 
(2£l.\  ~  „ 


f*  pressure  drop  across  the  cooler  (air  side' 
diameter  ratio  for  cooler 


Tlr  -  enthalpy  difference  across  the  interior 
stream  of  cooler 


v 

'■f  3, 
•  *  *  *  * 


sp.  vol.  of  working  fluid  at  WH3  exit 
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'a?\ 
— //j-HB 

AE\ 

,  P  X  ,'TT. 


V., 


fS  pressure  drop  across  the  ;iKB  ;  air  side' 

ft  pressure  drop  across  the  heater 
(inside  tubes) 


poiytropic  cooling  water  pump  e: 
sp.  vol.  of  cooling  water 


.c  ler.cy 


,  a?' 


/5k' 


Values  for  \  ?  /’Glr  and^/Clr  are  estimated  from  the 
following  relations: 


\  o  \ 

TV  C,,- 


■  J*P\ 


P  ,  - 


\ p  ,  jp  '  A?  . 

D  /  .«  •-‘5  x  ®  /  ^3*1  ■  °  ' 


~  /'ai"^T.  .  0"7  i 


A  n  /c.r 

Subscripts  1  and  2  refer  to  quantities  inside  and  outside 
the  tubes  respectively.  A„-,_  is  computed  as  follows: 

w  XX 


(22. 

I  O/r., 


,  .  4t 


;  -  =< 


—  •  \  __  /  .  ,  !  ! 

Aci-’a'V  '  =  ~ 

where  r  ,  ?  are  as  previously  defined  and 


-A  \ 

-  b  />;- 


—  /r 


.8 


\  /  1  /■ 


'  •  ••»  /  , 


( ,J  ' 
{ 


w  \ 

A  Vi  /r^/ 

'I*'/ ;/■. -  \  /  =* 


-  /A*, .  N  •  c  d  ^  o  yV  ■  ^  .  / 
'  \  /’  v  . 


y  V  /  •  / 


<  -  *  *  'i  /Uv  *  \  ’ i 


.  •  \w  /  ■ 

r  /  ’  )  1 

\  /  ''"i/ 


.  r'/h  ] 

\  /  'ri/ 
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he  volume  is 


Volcir 


,  /  =■ 

Acir  (.  p  /ci: 


d.  Heater.  The  heater  is  assumed  to  he  a  shell  and 
tube  heat  exchanger  with  the  wcri-cing  gas  in  the  cubes  and 
combustion  gas  out. 


'he  follow ins  are  defined: 


V  ^  w-. 


\  r  ,  ^  — - 

«  •  I  <4  <«4 


diameter  ratio  for  heater 

enthaphy  difference  across  the  i: 
stream  of  heater 


Values  off  A?'  and  (Th) 

\  /n  —  i  •  ?  /  n~2  *  -  u  • 

estimated  from  the  following  relations: 


\  r  /  <vr-  » 


:<c 


a?  \ 


\  *  /  r 


\  P  J h Ti-J  ■  \  2  /  ntr 


where 


^  /A  \~1,75  ( a  \"  f  -  \  p  /  ' 
^  (  /'wJ  1  /Ay  I'V/vi  /?4/ 


(37, rr  s  .04?TC  +.2/+ 


Subscripts  1  and  2  refer  to  quantities  inside  and  outside 
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Boiler  volume  is  estimated  using  heat  transfer  relations 
and  volumes  from  similar  U.S.  Navy  CCGA3  (Combined  las  Turbine 
ana  3 team)  designs  (ref.  ii). 


The  following  are  defined; 


ave.  specific  heat  on  steam  side 


ave.  specific  heat  on  air  side 


Number  Transfer  Units 


r.umoer  ci  oasses 


Heat  transfer  area 


overall  heat  transfer  coefficient 
total  volume 


From  a  DT.VSHi CC  (David  Taylor  Naval  Ship  Research  and 
Development  Center)  study  (ref.il)  the  following  relations 
are  obtained; 


CiV  f  ,V 

l.N/tX 


/n  ^  up 


c-r'« 


-  /'o.  ,,  Ci?  "  /  \  ,y.t  „  | ' 


N  U  =  (*  J  < 1  [i  -  CC6  -  £ 


'/  ->*s 


AI7U  *  C.’d/N 


For  similar  designs  assume  Av»,  =<  V„ 


r? 


V-  CMIN-  NTU, 


'll 


cmint 


II 


vT  ^  /v— 


•  I  —  Utt 

lx 


cm  IN-  ::tut 


'MTV  V*!T  , 

\uiiUltfT  .  <  -  w  —  —  / 

'  -k  X  a.  J.  / 


The  iera;V., 


an  be  evaluated  from  known  data 


'  '  Tl  7  T 

J  — —  ' 
X-  -  X, 


on  crevious  designs  to  rave 


V.  .... 

.JftJ 


=  11.53  *  CMIN  *  NT” 


But  it  should  ce  noted  that  the  above  equation  assumes 

AC  \ 


and  boiler  cower.  If  as  oeror 

/  .« r.3 


constant 
then  from  heat  transfer 


s  hefars  7n'_  “v. 


*V  n 


,4^-  *  >  «  Aos 


H  ♦ 


:hus 


■/  °<  r^.  f  An  $ 

Correcting  the  above  equation  for  power  gives 


V,/K3  =  .001?  CMIN  *  MTU  * 


\  - 


-  4./ 

AP\ 


As  stated  previously  Vol  °<  l~p~y 
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Correcting  again 


pressure  drop  yields 


- 


.0CC45  CI-TIH  *  NT’J  -‘i'A,  *■  >.J  ~~ 


o 


f .  Condenser .  Che  condenser  is  assumed  to  be  a  shell 
and  tube  crossflow  heat  exchanger  with  salt  water  inside  the 
tubes  and  steam  outside. 


Che 


■\r 

.  t 


A\ 


v 

c 


following  are  defined. 

inside  tube  diameter 
number  of  tubes 
length  of  tubes 

C-ravitatior.al  acceleration  constant 
pressure  drop  on  the  salt  water  side 
friction  factor 
density  of  cooling  water 
velocity  of  coding  water 


A„  -  cross  sect,  free  flow  area  on  water  side 
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w 

r * 

flew  rate  of  cooling  water 

q 

heat  flew  in  condenser 

w 

•n  x;  specific  heat  of  cooj.rrg  wa -e 
c  '  * 

t> 

o 

tamo,  difference  across  the  ceolm 

X 

steam  quality  entering  condenser 

i 

-  o 

er.thapy  difference  between  sat.  li 
sat.  steam 

A, 

mass  flow  of  steam  entering  ccncsr 

From  'oasic 

heat  transfer  relationships 

A,_jm  =  l! 

where 

A  \  . 

N)  “*  -W'  - 

.  »  ^  -s 

therefore , 

■)  1  5  \ 

:v.  DESIGN  C?  A  2  5,000  H?  PLANT 


A  proposed  design  of  a  25.CCC  H?  coal -fired  conbir.ee 
cycle  was  completed.  Goal  combustion  is  accomplished  using 
a  cressurized  fluidized  bed  (PFB)  .  Design  parameters  were 
eventually  chosen  after  investigating  the  effect  of  changing 
pressure  ratios,  varying  “T"  ,  -  ir.  the  crayton  cyc_e,  varying 
£  ; regenerator  effectiveness),  different  ccmcresscr  inlet  te 

(1.  )  and  steam  condensing  temp.  and  variations  in  Pray 

wonting  pressure  (.<?)  and  saturation  pressure  ). 


a.  The  Fluidized  Bed.  The  oressurized  fluidized  bed 


operates  at  a  pressure  of  5  12  atmospheres. 


“  c.  r  t:  '  * o  '  -r 

—  —  w  vU.  —  — 


maintained  by  an  exhaust  driver,  compressor.  For  proper  com¬ 
bustion  incoming  air  must  be  preheated  by  the  recuperator  to 
at  least  70C°7.  The  bed  material  is  primarily  limestone  whic 


serves 


ipture  sulfer  in  the  coal  and  thus  minimize  SC-, 


emissions.  A  bed  temperature  of  1?0C°F  is  generally  used  to 
heep  NO  emissions  at  an  acceptable  level. 

The  fluid  bed  itself  is  about  1  ft.  deep.  Preheated 
combustion  air  enters  the  bed  via  a  perforated  bed  plate  ’or.de 
r.eath  the  bed.  Pressure  dreo  across  the  bed  is  assumed  to  be 


1C5. 


The  PF3  is  used  to  heat  the  working  fluid  in  the  closed 


3rayton  cycle.  The  working  fluid  leaves  the  regenerator  and 
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Assumptions 


The  computer  program  was  run  with  one 


following  ir.pu 

t  va 

lues : 

0.  !?' 
1-'"' 

= 

.33 

■>;.  fp'. 
i  *  *  - 

= 

•  91 

-- 

= 

.33 

^<slSJ 

= 

.30 

fis 

1  • 

= 

.33 

■AL\ 

= 

• « c 

,"N 

/  ...-3 

= 

•  9^ 

-n 

•  < 

1 

= 

5^0.0  °R 

•  1, 

= 

2360.0  °? 

*  «  f7 

= 

51-0.0  CR 

/A : 

= 

-35.0  °R 

> «  i" 

= 

x^5*^  to* 

•*  7?  ^  <T> 

= 

200. 0  psi 

(want  at  least  =C  psi  greater 
than  W?) 

Power 

= 

250CO  H? 

T» 

~  lh 

- 

535.0  °R 

m 

*7h 

?6o  °R  (limited  by  acid  corrosion  in  the 
turbine ) 

= 

.35 

s 

.39 

- 

.90  (11* 

excess  air) 

In  addition  i 

t  was  assumed  that  scoop  injection  could  be 

used  for  both 

the 

cooler  and  c 

cndenser  at  speeds  greater  than 

10  knots,  ther 

efore  powor  requirements  for  salt  water  circuia 

pumps  were  excluded  ir.  the  calculations. 

c .  Results .  The  computer  output  is  provided  in  Appendix 
3.  Note  that  the  saturation  pressure  for  pressure  ratios  of 
2  and  3  were  below  the  2CC  psia  specif ication .  This  is  due 
to  the  fact  that  the  low  pressure  stream  leaving  the  regener¬ 
ator  is  at  a  temperature  below  the  saturation  temp,  at  ICCpsia. 

A  oiot  of  C-  vs.  r  and  Vol*  *  ,  vs.  7T„  is  shown  in  figure 

c  total  c 

a 


V  a  ^  **  V\  Q 

,  1  O  wCu  «  wtiC 

volume  increases 

sharply  above 

-  - 

due  largely  to  the 

rapid  growth 

in 

the  waste  heat  tcile 

V* 

Maximum  occurs 

at  around  .  1 

6.  To  investigate 

the 

effect  of  [^rjm  or. 

V-  and  Vclt, 

the 

program  was  run  for 

a-o  - 

— — .  =  jU  06 

.06,  .10,  and 

.1 

2  between  TT,  of  ** 

and 

The  results  plotted  or.  Figure  9  shew  that  significant 
gains  in  y_  can  be  realised  with  very  little  increase  in  7„ 
for  the  range  of  A- -6  pressure  ratio.  The  effect  of  changing 
Z-,  (regenerator  effectiveness)  in  this  range  with  a  - 

.05  was  further  investigated. 

The  results  listed  in  Table  1  show  no  increase  in  effi¬ 
ciency  as  Z^  is  reduced.  Furthermore,  volume  increases  sharply 
as  a  greater  percentage  of  the  work  is  being  done  by  the  Sank in e 
cycle  as  Z^  decreases. 

Thus  for  greatest  efficiency  with  minimum  volume,  cycle  operation 
should  be  at  the  highest  Zp  practicable. 


I  r}  'O 


'To  investigate  the  effect  of  changes  in  T,  (compressor 
inlet  ten?,}  runs  were  made  or.  4  different  temperatures  and 
the  results  shewn  in  Figure  10.  As  expected  efficiency  de¬ 
creased  with  increasing  T1  with  very  little  change  in  the  range 
of  interest  W.j  Volume  or.  the  other  hand  decreased 

slightly  for  less  than  3-5  'out  increased  significantly 
above  that.  Thus  for  best  efficiency  ar.d  least  volume  T. 
should  be  beat  as  close  to  cooling  water  tent erasure  as  tcssibl 


Finally,  changes  in  working  pressure  of  the  Bray  ter.  cycle 
were  investigated.  Tecreases  in  wcrbir.g  pressure  may  require 
additional  volume  of  3raytor.  cycle  components  yet  at  the  same 
time  allow  a  decrease  in  ?. aricine  cycle  steam  pressure  (recall 
pressure  differential  between  boiler  streams).  This  will 
in  turn  increase  cycle  efficiency  yet  further  increase  volume. 


These  results  are  corn  out  in  Figure  11. 


d.  System  Selection.  In  accordance 
criteria  the  highest  thermal  efficiency  c 
(Figure  11).  Thus  it  is  desired  to  selec 
imum  volume  at  </r  =  -432. 


with  the  design 
or.siaered  is  0.437 
t  the  plant  with  rni.n- 


Figure  11  shows  that  only  those  alter r. 
cycle  working  pressures  of  7Cpsi  and  lOCpsi 
requirement.  Gf  these  two  the  ICOpsi  case 
at  desired  efficiency. 

However,  the  problem  does  not  end  here 


atives  with  Bray  tor. 

meet  the  efficiency 
has  minimum  volume 

.  If  one  locks  at 
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On 


the  steam  quality  ieavir.g  the  turbine  at  the  77^  =  5,  = 

100  psi  point  it  is  of  the  order  of  0.21.  This  is  a  little 
lower  than  normally  desired  for  design.  In  fact,  the  steam 
quality  of  all  alternatives  thus  far  considered  within  the 
range  of  interest  is  low.  If  a  steam  quality  of  order  .3?  i; 
desired  the  current  design  ?t.  oan  change  in  cr.e  of  3  '•"'ays 
as  follows: 

1)  Raise  steam  condensing  temo 

2}  Lower  minimum  saturation  pressure  (?„.„) 

3)  or  a  combination  of  the  above  two 


Through  trial  and  error  it  was  found  that  case  1  gave  the 
greatest  benefit  with  the  least  penalty.  Computer  runs  were 
made  in  the  1-5  77^  range  with  the  steam  condensing  temp 
raised  to  560°R  for  both  the  70  and  ICCpsi  alternatives.  The 
results  shew  a  substantial  increase  in  quality,  slight  decrease 
in  efficiency  and  very  little  cnar.ge  in  total  volume  (Table  2)  . 


(psi)  __ 

Wording  Pressure  ■'  c 

u 

X 

Vo1-.  a-.h 

70 

^•5 

.180 

.861 

913.3 

7C 

5.0 

.182 

.871 

mo.  6 

100 

5.0 

.179 

.362 

956.5 

Table 

iu. 

The  70 

psi  designs  have  better 

steam 

quality  and  0  *  than 

/  w 

the  best  100 

psi  case.  In  addition, 

-  i± 

5  T„  the  volume  is 

lower.  3e  tween  the  two  70  psi 

cases 

above 

a  trade-off  exists 

between  9  t 

and  V 0 1 * .  Assume 

w 

that  volume 

comparable  to  the  case 
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os 


ICO  osi 


i  -  > 


"c  ~  •'* 

is  desired. 


,  =  5i-o°R  discussed  accve  ; Figure 


If  this  is  the  case  then  the  new  design  zcir.t 


is  at  a  Brayton  cycle  working  pressure  of  70  psi, 


=  5oC°H.  The  results  of  this  design  are  sumariced 


x.ne 


1 


following  cor.Duzer  or  in  tout : 


"5t5  -  vrnr^v  s*'?  ‘ 
^rv.A-  - 


- —  .  fl-r-.  _ 

.-.niwCii-r.-  --  ---  -  -  - 


WORKING  PHIS SURE  =  70.0  ^ 
OCTAL  PRESSURE  IRC?  =  •  d? 
MASS  PLOW  (  L3M/ S  ,  -  2  C  o  .  9  3 


‘T'TT 

MP( IEG 

3  > 

•  «.  / 

r*,n  *  ■— 

w  v  i* 

1 

5^0.0 

2 

376.0 

3 

1W03.9 

jL 

2CoO .  0 

** 

1-75-9 

* 

0 

9-3 .0 

?P 

a  ^6 . 2 

7 

309-2 

"'OS 

m  -m 

2  HU 

•rr.r  '  rv  ' 

*  ^  ^  ^  ^  V  4 

r'  *  r\ 

•  w  L  ^ 

.95c 

TUR3INE 

5  •  5 

SOCLE?. 

.009 

.^90 

4  3  *5  ' 
-7  J  •  ~ 

COMPRESSOR 

.  GCCo 

29.7 

■  *TT  '  ■nT  E> 

3CILER 

.030 

PRESSURE  RAT 

IC=w. 70 

,  .iVT’l  "*  ^  ~  ~  T  T  r  T  ,  -L  C  ' 

spc(L3M/H?  -kr)=c 

l  'w  _  w,  iVA-j  V  -  *  *  -  i 

POWER (HP ;=25CCC. 

work  era:  by  si:-: 


’.5  i??v?=c.9l 


29.7  Z??(?)=c.3: 

EE?  =  0  ,9«- 


n,’CmAL  PRESSURE  IRC?  =  .-3 
AIR  PLOW ( L3M/5 ) =  36.35 
FUEL  PLOW \ L3M/S ) =  3- Co 

STATE  TEMP C DEG  R) 


IK 

535-0 

335.2 

3K 

iLu-O.  5 

W-K 

U-SC9 .3 

5H 

1503.9 

6  H 

9^5.0 

7H 

760.0 

C  MFC NENT 

IP? 

IKD 

CO 

RECUP 

•  313 

.650 

1.9 

TUR3INE 

0 . 7 

HEATER 

.019 

.L30 

31 .* 

OPPRESSOR 

10 .  s 

RANKINE  CYCLE  j  SATURATION  PRESSURE  =  120.0 
STEAM  PLOW ( L3M/S )  =  5.69 
PLOW  RATIO  =0.0275 

STATE  TEI»t?(LEG  R)  PRESSURE  (  PS  I ) 

IS  560.0  C.9 

2S  575-0  dO.O 

S?  301.2  120.0 

**S  913.0  120.0  _ 

is  56O.O  0.9  LuAi.=0.c6^ 

CO  M?C  NENT  7CL( CU _ ?T ) 

TURBINE  6.5  EPP( IS) =0.33 

SON TENSER  13u-2 


Ct< 


»  ^  Or  **  *  m  • 


The  application  of  this  system  to  Mavy  ships  has  nary 
implications .  Cr.e  of  which  is  the  fact  that  neither  Blosed 
Bray  tor.  ncr  combined  cycles  have  been  previously  used.  However, 
there  is  a  considerable  effort  at  this  time  to  get  CCGA3 
(Combined  C-as  Turbine  and  'Steam)  systems  cn  board  Havy  ships . 

\  i  ■■pH  “tlfh  Vi ^  ^  rv  '  T  ■“*'  *  o  o  _  "n  r>  0  3^035'^ 

Bray  ter.  combined  cycle  when  equipped  with  a  multi-fuels  heater 
offers  the  additional  flexibility  of  burning  the  cheapest  cr 
most  readily  available  fuel. 


These  systems  are  all  limited  to  burning  oil  fuels 


Advantages  of  the  Closed  Braytor.  combined  cycle  are: 

1)  The  Braytor.  cycle  operates  at  temperatures  substan¬ 
tially  below  the  material  limit  than  ir.  the  open  cycle  gas 
turbines  resulting  ir.  cheaper,  mere  reliable,  and  longer  life 


components . 

2)  low  pressure  steam  systems  are  well  proven,  ar.d  have 
significantly  less  water  chemistry  problems  than  higher  pressure 
ones . 

3)  As  lead  decreases  and  fuel  flow  is  cut  back,  the  con¬ 
trol  valves  (Figure  12)  between  the  compressor  outlet  ar.d  cooler 
inlet  car.  be  opened  as  necessary  (keeping  max.  temp  and  speed 
constant)  so  as  to  ensure  that  the  vector  diagrams  of  all  sec¬ 
tions  of  the  turbomachinery  blading  remain  unchanged.  Turbine 
ar.d  compressor  efficiencies  will  remain  essentially  the  same 
producing  high  <}  *  at  part  load.  This  requires  use  of  a  CR? 
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(controllable  reversible  pitch)  propeller  to  maintain  constant 
shaft  speed  at  part  loads. 

Disadvantages  of  the  Closed  Braytor.  combined  cycle  are: 

1)  Higher  overall  system  weight  and  volume  due  to  the 
additional  heat  exchangers  and  piping  required  for  the  Closed 
Bray ton  cycle . 

2)  A  complex,  integrated  control  system  is  required  to 
coordinate  control  of  the  combined  cycle,  fuel  flow,  and  C?.? 
at  part  loads. 

3)  Fluidised  bed  combustion  is  not  a  fully  developed 
engineering  concept.  Land  based  prototypes  have  been  built 
and  operated  with  limited  but  increasing  levels  of  success.  A 
shipboard  proto  type  is  perhaps  5  years  in  the  future  with  pro¬ 
duction  of  a  mar inis ed  version  at  least  10  years  away. 

k)  If  coal  is  to  be  used  as  the  primary  fuel,  it  should 
be  nomad  that  it  has  a  lower  heating  value  per  ton  than  oil. 

The  result  is  a  larger  potentially  more  expensive  (greater 
acquisition  cost)  ship  for  the  same  military  mission. 

5)  Reaction  time  with  coal  fired  fluidised  bed  combustors 
to  load  changes  are  generally  slower  than  that  for  oil  heaters 
necessitating  in  a  separate  fuel  oil  system  for  rapid  increases 
in  load. 
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VI.  CONCLUSIONS  AMD  P.ZCGMMEMDATICMS 


The  Closed  Brayton  combined  cycle  offers  an  attractive 
alternative  system  for  many  shipboard  applications.  As  with 
the  first  Arab  oil  embargo,  subsequent  cut-offs  of  Arab  oil 
could  result  in  decreased  steaming  days  and  a  lower  overall 
readiness  of  the  U.S.  Navy's  defense  posture. 
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Paramount  in  the  engineering  of  a  successful  multi-fuels 
propulsion  plant  is  the  development  of  a  mar ini ced  fluid  bed 
combustor.  Che  design  of  heat  exchangers  and  other  turbo- 
machinery  used  in  the  combined  cycle  are  well  within  the  state 
of  the  art  of  present  technology. 

It  is  thus  that  the  author  concludes  the  following: 

1)  The  design  of  a  modern  coal  burning  combined  cycle  is 
within  the  reach  of  present  technology. 

2)  This  combined  cycle  would  offer  high  efficiency  at 
both  rated  load  and  part  load  due  to  the  unique  controlling 
aspect  of  the  Closed  Brayton  cycle. 

3)  '/.'eight  and  volume  of  a  Closed  Brayton  combined  cycle 
would  be  greater  than  that  of  an  equivalently  rated  CCGA3  type 
system. 

4)  Coal  offers  the  benefit  of  a  readily  available, 
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relatively  cheap  fuel  at  the  expense  of  increased  ship  size 
and  acquisition  cost  due  to  lower  2fU  heat  value  per  ter.  when 
compared  to  oil. 

The  author  maices  the  following  recommendations : 

1)  Development  of  a  land  based  test  site  of  a  Closed 
Brayton  combined  cycle  equipped  with  a  fluid  bed  heater  sr.ouid 
be  considered. 

2)  A  pressurized  fluidized  bed  combustor  suitable  to 
marine  applications  should  be  developed. 

3)  A  shipboard  prototype  system  should  be  integrated 
on  a  U.3.  Navy  ship  to  attempt  to  quantitatively  define  the 
shipboard  impact  of  the  system. 
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APPENDIX  A 


The  computer  program  performs  a  thermo dynamic  analysis 
and  component  volume  computation  for  various  specified 
Brayton  cycle  compressor  pressure  ratios. 

The  program  accepts  as  input  the  following  quantities: 


%  ( ? ) 


/-  (?) 


•  v 

/  h  P  \ 

[~J7 

/  IK 3 


polytropic  compressor  efficiency  (Brayton 
cycle} 

polytropic  turbine  efficiency  (Brayton 
cycle ) 

regenerator  effectiveness 


isentropic  pump  efficiency  ( Panic  in  e 
cycle) 

total  pressure  drop  around  Brayton  cycle 
loop 

efficiency  of  the  waste  heat  boiler 
compressor  inlet  temp.  (Brayton  cycle) 


■u 


bur  bine  inlet  temp  (Brayton  cycle' 


*13 


A 


Power 

m 

~lh 


/«"( is) 


steam  condensing  temp. 

min.  temp  diff.  between  streams  in  V.'KB 
(at  ??  and  Tg) 

min.  saturation  pressure  in  the  Rankins  cycle 

net  power  out  of  the  combined  cycle 
compressor  inlet  temp.  (heater  cycle) 

turbine  exit  temp,  (heater  cycle) 

polytropic  compressor  efficiency  (heater 
cycle) 

isentropic  turbine  efficiency  (Rankire 
cycle) 
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\ 


a 


w? 


poiytropic  turbine  efficiency  (heater  cycle; 
equivalence  ratic 

working  pressure  in  the  Brayton  cycle. 


Output  is  as  shown  in  Appendix  3. 

The  main  program  performs  the  thermodynamic  analysis  cf 
the  closed  Brayton  cycle  and  the  ParJcir.e  cycle.  Functions  per¬ 
formed  by  the  various  subf unc tiers  are  as  follows: 

1)  Subroutine  HZ A TIP  performs  a  thermodynamic  analysis  of 
the  heater  cycle. 

2)  Subroutine  BHD  computes-—1  for  the  regenerator,  re- 


cuperatcr,  cooler  and  heater 


±p 


3)  Subroutine  DPP  calculates  ns- for  the  regenerator,  re 


cuperator,  cooler, 

4)  Subroutine 

5)  Subroutine 
5)  Subroutine 


*?  > 

f  / 

3  x 

^  J 

9) 

10) 


1 


1) 


4 

4. 


2) 


tables . 


Subroutine 
Subroutine 
5  ubroutir.e 
Subroutine 
3  ubroutir.e 
Subroutine 
Reproduced 


heater  and  boiler. 

TUP7CL  finds  turbine  volumes. 

CCMPVCL  finds  compressor  voi-ur.es. 

PZGZN  computes  the  regenerator  volume. 
PZCB'P  computes  the  recuperator  volume. 
CCCLEP  computes  the  cooler  volume. 

CO M3 US TCP  computes  the  heater  volume. 
WK3VCL  computes  the  waste  heat  boiler  vo 
CCNDVCL  computes  the  condenser  volume. 
STMT A 3  is  a  computerized  version  of  the 
from  ref.  15* 


e  . 


s  team 
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DIMENSION  C R ( 3 0 )  ,SATLIC(3 )  ,STEAM(3 }  ,7??C2)  ,WR(11) 

REAL  MCP ,NT? .HP I S ,NTI S , MWHB ,N7,Ll  ,L2  ,KO 
R  E  A  L  M  C?  H ,  N  T  P  H ,  M  8 ,  M  S ,  M  A ,  M  F ,  N  H , I ,  I T ,  .<  1 1 ,  M  P  P 
REAL  LHV,M,M1,M2,M7U 

COMMON  /STM/PRESS, TEMP .QUAL.SATLIQ, STEAM 

COMMON  /HTR/NH,PRH,DPPH,T2H,T3H,7AH,T5H,~oH,;!A,MF 

COMMON  ,/D IA/3HDREG.3HDCL.7 ,3 HOHTR ,3HGP.EC 

COMMON  /PRES/0PPyH3,DPPREGfDPPCLR,DPPREC,3PHl,3PH2,3H0H7 

COMMON  /CCMB/QPPRG,3PPRC,3PPHE,2HDHR,3PME2 

READ  ( 5 , 1 0 )  MCP , NT? , ER , NP I S  ,NT 1 3  ,D ? ,  NUH3  , 71 , 74 , 71S  ,3  7 ,  ?M I  n , :< 

FORMAT  (7F4.3, 5F7 .1,12) 

READ ( 5 , 1 5 }  PUR , 71 H , 77 H , M C? H , NTP H ,? H I , „P 
F0RMA7  ( F3 . 1 , 2F6 . 1 , 3F5 . 3 , FS . 1 ) 

REAO(5,20)  (CRH), 

FORMAT  ( <K>F5  .2) 

R3  .063544 
OC  200  >l,i< 

PR=CR(  I ) 

Z3C?  1(  T1  )•*■(  R/NC? )  ■*LOG(?R ) 

T2aEXP(Z/ .2475) 

CALL  ITER  (72 ,Z) 

7R3!/  (1-OP)  /'PR 

Z=CP 1  (  74 ) *R*N7?*LCG( 7R ) 

7S3EX?(Z/.2475) 

CALL  ITER  (T5,Z) 

IF(T2  .GT.  T5)  GO  TO  30 
763T5-ER*(T5-72) 

T33T2+E?*( 75-72) 

GO  70  40 
76  3  75 
73*T2 
LIaTl 
L23I 122,0 

IF(L2  .GT.  To)  GO  TO  50 
GO  TO  60 

L2aT6-(T6-Tl) /10 . 

TPP(l)aL2 

WR1»0. 

Ca(L2-Ll) /10. 

DO  70  >1,11 
TSFaTPP(J)-OT 
T3S3T5-0T 
T«TSF-4oO 

CALL  STMTA8  (2,7,0.) 

PSF»PR£SS 

HSFaSATLIQ(2) 

TaT3S-460 

CALL  STMTA8  (7.PSF.7) 

H3SaSTEAM(2 ) 

S3S-STEAM(3) 

T*TlS-460 


CALL  STMT A3  (2,7,0.) 

?4S=PRE$3 
V 1S=SATLIQ ( 1 ) 

H1$=SA7LIQ(2) 

SIS— SATL IQ ( 3 ) 

CALL  STMTAB  (9,?4S,S3S) 

H4SIS»STEAM(2) 

H4S»H3S-MTIS*(H3S-H4SIS) 

CALL  STMTA3  (8,?4S,H4S) 

X=QUAL 

Y4S*STEAM(1) 

H2SIS»HlS+VlS*(PSr-P4S) 

H2S*HlS+( H2SIS-H1S) /MP I S 
T2S=T1S+15. 

B-NWHB* ( C?2 ( T6 ) -C? 2 ( TPP { J ) ) ) / ( H3S-HSF  ) 

Z«AC? ( TPP( J ) ) *T?P ( J ) -3’ ( H  SF-H2S ) /HVH  3 
T7  =T?P ( J ) 

CALL  I7EP2  (T7,Z) 

WR  ( J )  s3'*(  H3S+H1S-H4S-H2S 1 
IFfJ  .EQ.  1)  GO  TO  58 
IF ( T7  .LT.  Tl)  GO  TO  75 
IF(PSF  -LT.  PNIM)  GO  TO  75 
IF(WR(J)  .GT.  WR1)  GO  TO  63 
IF(X  .GT.  .37)  GO  TO  75 
53  wRl=WR( J ) 

TPP(J+1) =T?P(0 )-C 
70  CONTINUE 

GO  TO  76 

75  IF(ABS(TPP(J)-T?P(u-l)  )  .LT.  .1)  GO  TO  76 
li-TPP(J) 

L2aTP?(u-I) 

GO  TO  50 

7  6  WB-C? 2( TA ) *C?  2( Tl ) -C?  2( T5 ) -C? 2( T2 ) 

QIM«CP2(T4)-CP2(T3) 

MXMB+MRU)  )/QIM 
MB».7Q7*PWR/ (WB+WR(J) ) 

MSaB*MB 

QB3MB*QIN 

CALL  HEATER  (QB  , T1H,T7H,MC?H, MTPH, PHI, T3, 741- 
CALL  OHD(PR) 

OHOHTsQHDHTR 

K»1 

77  CALL  0PP(T6 ,T7 ,?R ,DP,£R,T2,T3 ,DP°H,74H ,T5H ,T4,?RH,UP ,MB ,MA 

*  ,mf,;<) 

Ga1.33 

N*G/(G-(G-1. )*NT?) 

PaWP*(PR-OPPREG) 

V3l./ARHQ(T4,P) 

VMNTURB=Q.Q 

CALL  TURBVOL  (N,MB ,V,TR,T4 .YMNTURB ) 

G-1.39 


62 


>J*G/'  ( G-(  3-i .  i  *N7?H ) 

?*14. 7* (PRH-OPPREC) 

7*1./ARHQ(T5H ,? ) 

M*MA+MF 

VHTRTUR8*0.0 

T?H*1./ (l.-OPPH).'PRH 

CALL  TUR3VOL  ( N ,  M ,  V ,  7R  H ,  75  H ,  VHTRT’d  R  3 ) 

7*T1S/T3S 

?*P4$,'PSr 

M*I. / (1 . -L3G(T) /L0G{ P ) ) 

YSTMTUR3-0.0 

CALL  TuRSVCL  ( N , MS , V 4S ,? , T3S , V STMTUR3 ) 

Gal  .4 

M*G/(G-(G-1.  ),'lC? ) 

7*1./ARH0(71,WP1 

7MNCQMP*0.Q 

CALL  COMPVCL  1  7, MS  ,7  ,?R .7MNC0MP ) 

7*1. / ARHO (535.  ,14.7} 

VHTRC0MP*0.0 

CALL  CCMP7CL  (U.MA.7  ,?RH,7H7RCCMP  ) 

0*  .125 
V REGEM* 0.0 

CALL  REGEM  (T2,T3 ,75 ,75 ,*P  ,3 ,M3 ,?R ,3?PR£S,0HDR£G, YRE3EM ) 
7REC'JP*0.0 

CALL  RECUP  !T2H  ,~3H ,75H ,7SH  ,?RH,3??H,3  ,MA  ,MF  (3?cREC  ,7R£CUP ) 

3*  .25 
YCLR-0.0 

CALL  COOLER  ( 77  ,T1 ,  h?  .O.MTP.MB  ,3H0CLR,0PPCL.R,7CLR) 

OPPRG*OPPREG 

OPPRC=OPPREC 

DPPHE=QPPH 

0H0HR-0HDH7R 

0PHE2-0PH2 

VH7R*0.0 

CALL  COMBUSTOR  ( 73 , 74 , 74H , 75 H , *P , P R ,3 , MB ,? R H , MA , MF ,7H7R) 
VyHB*0.0 

CALL  WHBVOL  (MB ,MS,76 ,77 ,MWHB ,H2S ,H3S ,3PPWH8 ,VWH8) 

VC0M0*O.0 

CALL  CONDYGL  (MS,X,7C0ND) 

V0LT-VMNTURB+VHTRTURB+VS7MTUR3+VMMC0MP+VHTRC0MP+VREGEM+VRECUP 
*  +YCLR+VHTR+YWHB+VCOMO 
IFOC  .EQ.  1)  GO  TO  78 
!F(VOLT  ,GT.  V0L1)  GO  TO  31 
IF ( (DPPCLR- .002)  .LT.  0.0}  GO  70  31 
78  K*K*1 

VOLl’VOLT 
V1*VCLR 
72*VWHB 
GO  TO  77 
31  70LT*V0L1 

VCLR*7I 
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/ 


7WH8=V2 


WRITE!5 ,32) 

PR ,  rtP  ,QP ,M8 

32 

FORMAT! ' 1 ' , 

'3RAY70N  CYCL 

*  PRESSURE3 ' ,F 

5.1/  IX, 'TOTAL 

*  (lbm/s)- '  ,F5 

.2//10X, 'STAT 

WRITE(6 ,33) 

33 

FORMAT (12X, 

'  1'  ,  5  X ,  F  6 . 1 ) 

WRITE(o ,34) 

T2 

34 

FORMAT (12X, 

'V  ,5X ,  F5 .1 ) 

WRITE(6 ,35) 

— j 
•  -/ 

35 

FORMAT (12X, 

'3' ,5X,F6.1) 

WRITE(6 ,36) 

T  4 
»  + 

36 

FORMAT (12X, 

'4'  ,5X,F5.r; 

WRITE! 6 , 37) 

~c 

i  ^ 

37 

FORMAT !12X, 

'  5 '  ,  5  X ,  F  5 . 1 ) 

WRITE(6 ,38) 

To 

00 

FORMAT (12X, 

'5' ,5X,ro .1 ) 

WRITE{5 ,39; 

T?=>( J  ) 

39 

FORMAT! 12 X, 

*?P ' ,5X,F6 .1) 

WRITE(6 ,90) 

T7 

90 

FORMAT! 12X, 

'7' ,5X,F6 .1// 

'<  '  Mi  '  ", 


WRITE (6 ,31)  OPPRE5,OHDREG,VREQEN1ER 

91  F0RMAT(7X, 'RS3EN'  ,4X,F4 .3 , IX, F4 .3 , IX, F7 .1 ,2X, 1 EFF= 1 
WRITE (6, 92)  VMNTUR3 ,MT? 

92  F0RMAT(5X, 1 7URSINE'  ,15X, F5 . 1 ,2X, 1 EFF(? )* ’F1 .2) 

WRITE (6, 93)  0PPCLR,DHDCLR,7CLR 

93  F0RMAT(6X, 'COOLER1 ,4X, F4 .3 ,1X, F4 .3 ,3X,FS .1 ) 

WRITE (6 ,94)  VMNCCMP.NC? 

94  FORMAT! 4X, 'COMPRESSOR'  ,14X,F5 .i,2X,  'EFr!?)»'  ,?4.25 

WRITE (6, 95)  DPH1 

95  FORMAT! 6 X,' HEATER1 ,4X,F5.4) 

WRITE (6, 96)  0PPWH8  ,NWH3 

96  FORMAT(6X,  ‘30ILER1  ,4X, F4 .3 ,15X,  'EFF«'  ,F4.2 /) 

WRITE (6  ,100)  PRH,DPPH,MA,MF 

100  FORMAT! IX,' HEATER  CYCLE:  PRESSURE  RATIO ' ,F5 .2/1X, 'TOTAL 

*  PRESSURE  DROP* ' .F3.2/1X, 'AIR  FLOW! 1 bm/s)= FS .2/1X, 1  FUEL 

*  FL0W(1 bm/s)a 1 ,F5 .2//10X, 'STATE  TEMP (deg  R)1) 

WRITE (6, 101)  T1H 

101  FORMAT! 12X, ' 1 h ' ,5X,F5 .1) 

WRITE (5, 102)  T2H 

102  FORMAT! 12X, ‘Eh1 ,5X,F6.1) 

WRITE! 6, 10 3)  T3H 

103  FORMAT! 12X, ' 3  h ' ,5X,F6.1) 

WRITE!6  ,104)  T4H 

104  FORMAT! 12X, *4h’ ,5X,F6.1) 

WRITE (6 ,105)  T5H 

105  FORMAT! 12X,'5h' ,5X,F5.1) 

WRITE (6, 106)  T6H 

106  FORMAT! 12X, ' 5  h 1 ,5X,F6.1) 

WRITE! 6, 10 7)  T7H 

107  F0RMAT(12X,  '7h'  ,5X, F6.1//5X,  'COMPONENT  OPP  3HD  70L(cuft 


6a. 


ERH*  i  T3H-72H )  /  (75H-T2H  'j 

WRITE (5  ,110)  3PPREC  .OHCREC  ,YRECU?,£RH 

FORMAT! ?X, ‘RECUR1  ,4X,  F4.3  ,IX,  F4 .3 ,1X,F7 .1 ,2X,  1  EFF=  ‘  ,F4.2) 

WRITE (6,111)  YHTRTURB  ,NTPH 

FORMAT ( 6 X,  'TURBINE'  ,15X, F5 .1 ,2X, ' EFF(?  }«  '  ,.-4 .2) 

WRITE (5  ,112)  QPH2 ,3H0HTR,7HTR 

FORMAT ! 5  X , 1  HEATER 1 , 4X , F4 . 3 , 1 X , F4 . 3 , 3  X , FS . 1 ) 

WRITE (6,113)  VHTRCOMP.NCPH 

F0RMAT(4X,  'COMPRESSOR1  ,UX,  "5 .1 , 2X,  'EFF{?  )* '  ,F<1.2/  I 
WRITE (6, 120)  ?SF ,MS,3 

F  0  RMAT  (IX,  1  RANK  I  ME  CYCLE:  SATURATION  PRESSURE*  1  ,F5 . 1,'IX, 

*  'STEAM  FLOWClbrn/s)*'  ,F5 .2/1X,  'FLOW  RATIO*'  ,-5.4/  '1CX,  'STATE 

*  TE>P(deg  R)  PRESSURE! psi) ' ) 

‘/WRITE ( 5  ,121)  T1S,?4S 

FORMAT! 12X, 'Is' ,5X,F5.I,3X,F5.1) 

WRITE (6  ,122)  T2S ,?SF 

FORMAT! 1 2X, '2s' ,5X, F5 .1 ,3X, F6 .1 ) 

WRITE (6 , 123)  TSF  ,?SF 

FORMAT! 12X, ' SF' ,5X,F6 .1.3X.F5 .1) 

WRITE (5 ,124)  T3S  ,?SF 

FORMAT! 12X, '3s' ,5X,F5 .1 ,3X,F6.1 ) 

WRITE (6  ,125)  T1S.PAS.X 

F0RMAT(12X,  '4s'  ,5X, FS  .1 ,9X,F6 .1 ,2X,  'Ql'AL*  '  ,F5 .3/  } 

WRITE (6,123)  V'W'riB 

FORMAT! 5 X,’ COMPONENT  VCL(cu  ft)  76X, 'BCILER' ,5X,F6 .1) 

WRITE (6, 129)  VSTMTUR3  ,NTIS 

F0RMAT(6X, 'TURBINE' ,3X,F3 .1 ,2X, 1 EFF ( IS)* 1 ,F4 .2) 

WRITE (6  ,130)  /COMO 

FORMAT! 5X, 'CONDENSER' ,4X,F5.1/) 

WRITE (6, 135)  NT 

FORMAT!  '1'  ,' SUMMARY: THERMAL  EFF*'  ,.-5.3) 

WRITE  15 ,136)  MH.PHI 

F0RMAT(9X,  "HEATER  EFF*',F5.3,'  (PHI* '  ,F4 .2, '  ) ' ) 

SFC*MF*3600 ,/PWR 
WRITE(6 ,137)  SFC ,YOLT,?WR 

FORMAT (9X,'SFC(lbm/HP -hr)*1 .F5.3/9X, 'VOL  TCTALtcu  ft)*1 ,F5 .1/ 
*9X, ’?0WER(HP)*’,F7.1) 

RANK  I N  E *W  R  ( J )  /W’B 
WRITE(6 ,138)  RANKINE 

FORMAT (9X,  "WORK  FRAC  3Y  STM  CYCLE* 1 ,F5 .3 ) 

CONTINUE 

*  STOP 
END 

FUNCTION  CP1(T) 

CPl*.2475*L0G(TM3.759E-5*TM2.553E-3*rt*2M4.i03E-12*T**3) 

RETURN 

END 

PJNCTION  C?2(T) 

C?2«.2475n-(1.379E-5*T-**2)+(1.702E-a*T**3}-(3.078E-i2*T**4) 

RETURN 
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END 

FUNCTION  AMT) 

TK=5T/9. 

AMU»TX**1 . 5*1 .  E- 6/  ; 7<+ 1 1 0 .  A ) 

RETURN 

END 

FUNCTION  ARHQvT,?; 

R®53.35 

ARHO-Pn-U./R/T 

RETURN 

END 

FUNCTION  AK'vT,?  • 

A= -532  5£- 5 
3=2A5. A 


T’<=5.'*T  9  . 

AK  >  ( ATT .  5 ) . '  ■  I .  *3*1 0 . ( - 1 .  M  '  ~< )  T< )  -A  1 3 . 

A'<  3  (  AX  0+  { 1 .  A6  A?  -  5*  ( 1 5 . 0  2 \ARHC  (  T  ,  ? ;  ;  1 . 2  2  'T . 

RETURN 
END 

FUNCTION  APR (T ,? ) 

C?=AC?(T) 

Xf1U»AWJ(T) 

XK«AK(T,P) 

APR«C?’*XM'J/XX 

RETURN 

END 

FUNCTION  AS  I  GMA  IT,?) 

A  S I GU  A=  1 7 . 2 1  A^I  5  2  2 2 . T  AX  I  "T ,  ?  )  )  T  - 1 .  A 1 )  T  .  A  c  T*U  i 
r  ( 1 . / APR ( T  ,?)*"*  .564)*  ( 1 .  /  (15  .C2*ARH0(  ~ ,?)  ;  **  .41 ) 
RETURN 
END 

FUNCTION  AC?(T) 

AC?3.  2475-3. 75?E-5*7»5.1C6E-3*T**2-1.231E-11*~~3 

RETURN 

ENO 
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SUBROUTINE  ITER  (T,X) 

10  T1*T 

7*EXP( :x+.24?5*L0G(T)-C?l(T) ’ /.2A75) 
IF( ABS(T-Tl)  .LT.  1.2-25  GO  TO  20 
GO  TO  10 
20  RETURN 
END 

SUBROUTINE  ITER2(T,Xi 
10  T1*T 

T*X/AC?( T) 

IF(A8S(T-T1)  .LT.  l.E-2)  GO  TO  20 
GO  TO  10 
20  RETURN 
END 
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SUBROUTINE  HEATER  ( OB ,TIH,T7H ,NC?H ,NT?H ,?H I ,T3 ,T4) 

REAL  L  HV ,  M  H ,  MF ,  MA  ,N  C?  H ,  N  7?  H ,  N 

COMMON  /HTR/NK ,PRH  ,DPP  h , T2H , T3H , T-JH , T5H ,  T6H ,MA ’F 

LHV*1290Q. 

r  A=  .0922 

H7H«AC? ( T7H ) *T7H 
H1H«ACP(T1H)*T1H 

N  H=  I  -  (  (1+FA~PH:)*H7H-H1H)/(FA*PHI*LHV) 

V,FSQ8/  (MH*LHV  ) 

MA*MF/(FA*PHI) 

T2H*0 .0 
T5H=0.Q 

10  C?T-AC?( {T5H+T"H)/2. ) 

CPC»ACP( (T1H+T2H  5/2. ) 

M=  ( i-f  A''?H  I )  ‘TC?”'rNT?H*  i  0?*',‘'7H-C?0'rTlH  ; 

DaC?C'*TlH*(C?C/NC?H-C?  TVJTPH) 

?RHa(N/D),,*(NC?H/C?C) 

TH2<lri*?RH*-*  ( CPC/NC?  H ) 

1+r  A*PHI}'*C?T*T7H+C?C'*Trrl'*  (RRH*’^  CPC/NCPH  )-I ) 

D*  i  I+F  A*PH  I )  *C?  T*T7H'*PRH*''  ( CPT/NTPH ) 

2P?Hal -(M/D )*■’'(  1 .  / ( C?T*NT?H)  ! 

TRH*(l-DPPH)irPRH 
TH6=T7H’r(TRH  )*•’'(  CPT'MTPH) 

IF( ABS(T2H-TH2)  .lT.  .1)  SC  TO  20 
IF(ABS( T6H-TH6)  .LT.  ,i)  00  TO  20 
T2H*TH2 
ToHaTH6 
GO  TO  10 
20  T2N-TH2 
T5NaTH6 
T5H3'iJt1QC. 

TAHaT5H+QB/(AC?(T5H/2+13C0. )’(MA^MF) ) 
T3H«T2H-f(l+MF,'MA)’r(T5H-T6H) 

RETURN 

EMO 

SUBROUTINE  DHD(PR) 

COMMON  /C I A/OHDREG ,2HDCL.R  .DHDHTR  .OHOREC 

0HDREG».557+.0933*PR 

IF (DHOREG  ,GT.  .95)  OHOREG*.95 

DHDCLR* .Q76+.Q92*PR 

DHDHTR3 .214+ ,048*PR 

DHDREC* .65 

RETURN 

END 

SUBROUTINE  DPP("6 ,T7 ,?R ,DP,ER ,T2 ,T3 ,DPPH,TH4,TH5 ,T4 ,PRH,P ,M,MA 
*.MF  ,,<) 

REAL  N ,MA ,MF ,M1 , M2 

COMMON  /PRES/CPPWH8 ,DPPREG,DP°CLR,DPPREC ,DPH1,DPH2,DHDHT 
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DTOT6-T7 

DP  PVJH83FL0AT  ( .< )  *  .002 
DPOR£G*( .135-.0294*PR)-0P/.l 

:f(oppreg  .lt.  .on  dppreg=.qi 

t.F(T2  .Eg.  T3 )  0pOR£G=0.0 
0PH2*(  .01Q8*PR-.0n*9P/.l 
II  *  ( T3+-T  A.  )■■*■.  5 
TT*(TH4+TH5)*.5 
Pl*P*(PR-0?) 

P2*14.7*{?RH-OPPH) 

M1=M 

M2=MA-t~MF 

XMUl*AMJ(71) 

Xt-iU2»AMU(TT) 

XR  HQ1=ARH0 ( ~i , P 1) 
XRH02«ARH0(TT,?2) 

3=  ( Ml  /M2  }*•*(- 1 . 7  5 1  *  (XM'J  1  /'(MU  2)  ~ 

3PHl=OPH2*OHDHT*,,3/3 

0  P  PC  L  R*0  P  • -0  P  PH  H  3-0  P  P  R  £  3-0  ?  HI 

0PPREC*QPPH-QPH2-  .1 

RE7URM 

END 


#  2  5 }  •*  'v 


hC’  XRH02'  ■* 


SUBROUTINE  TUR3V0L  (G,M,7 ,TR,T,VOL) 

REAL  M ,  M  ,.<  1 1 , 1 T 

SF=4.0 

111=850. 

PHIT=Q.6 
7AS0 .9 
CL=2.Q 
P>3. 14159 
RHGM=15,12 
A=9.S3Eo 
9--95.14 
AL?HA=1 .2 
Al-2.5 
Xl=  .307 
K Il=3.3SE-3 
THETA*0.3 

RT»SQRT(M*V/(PHIT*UT*P  IM1  .-YA**2) } ) 

V8=SQRT  ( 1 .  -TR**  ( I .  ,'G )  *  a . -VA**2 5 ) 

Cl«0.2* (Al+1.  ) *1 .  /V*(l . -VA**2 )  *1JT*‘*2*RT*'*2'*X  1/XI1 
C2a2 .* ( A-3*T ) / ( SF*ALPH A ) 
C3-2.*8*T*(1.-VA»*2)*»(G-1.  )/{SF*ALPHA) 
C4=-l.*THETA*RH0M*llT**2/2. 

H= (VA-VB ) /20. 

SUM*0.0 

X1*0.0 

Y1*VA 

00  100  J»i ,21 

M*C1*(  (l.-YD/d.+YDJ-PHIT^+Yl*^ 


i p  i  o  2 
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0*C2+C3*  ( 1 .  -Y  1**2)**  ( i .  -GDC4*  ( 1 .  -Y  1**2) 

DN=C1*(2 . *Y 1* (1 . -Y 1) / (1 . -Y D-(?HIT**2+Y 1**2 ) / (1 . *Y 1 ) -( 1 . -v 1 ; 
***(-2)*(l.-Yl)*(PHIT**2+Yl**2) ; 

0D=2.*Y1*(1.-Y 1**2 )**(-!. *G) *(3-1. )*C3-2.*Y1*C4 
X2*  .5* ( 1 ./ (M/D /** .5)*( (DN*0-0D*N)/0**2) 

IF (XI  .EQ.  0.0)  30  T0  50 
S0M=SUM+(X1»X2) , 2 . 

50  X1*X2 

Yl-Yl-H 
100  CONTINUE 
I7=H*SUM 

'/0L*2.*PI*R***2*CL*I7 

RETURN 

EMD 

SUBROUTINE  CCMPVOL  (G,M ,7 ,?R ,7CL) 

REAL  M ,  M  1 1 , 1 7 
Sr=4.0 
U 7*850. 

PHIT=0 .5 

VA*Q .6 

CL*2.0 

P 1*3.14159 

RH0M*15. 12 

ALPHA-1.2 

A1-12.Q 

XI*  .05 

■<U*4.93E-5 

7HE7A*0.3 

SIGMAY-12.54E5 

RT*SQRT(M*V/ (PHIT*U7*P I* ( 1. -VA**2) ) ) 

YB-SQRKI .  -?R**(-1./G  )*(1 .  -VA**2 ) ) 

Cl* .  2* ( Al  +  I . )  *1 .  /V*  ( 1 .  - VA**2 )  *117*  *2*R7**2*PH  1 7**2*X  l/KIl 
C2*-l*THETA*RHCM*U7**2/2. 

C3=2.*SIGMAY/( SF*ALPHA ) 

H*(VB-VA)/20. 

SUM-Q.O 

X1*0.0 

Y1*VA 

00  100  J*l,21 

N*C1*(  (1.-Y1)  /(1.+Y1)  )*(PHIT**2+.49*Y1**2)  / (P H 1 7**2 *.7 225* Y  1**2 
*) 

0*C3+C2* ( 1 . -Y 1**2 ) 

0N=C1*  (-1.  *(PHIT**2  +  .49*Y1**2 ) / (P H I 7**2 +.7 225* Y 1**2 ) *(1 .+Y1 ) -( 1 
*. -Yl)* ( 1.+Y1)** (-2)*(PHI7**2+ .49*Y 1**2) / (PHIT**2+ .7225*Y 1**2)+ 
*.98*Y1* (1 . -Yl) / (1 .+Y 1) / (PH1T**2  +  .7225*Y1**2 )-l .445*Y 1* (1 . -Y 1)  / 
*(1.+Y1)*(PHIT**2+.49*Y1**2)*(PHIT**2+ .7225*Y1**2)** (-2) ) 

00--2. *C2*Y1 

X2».5*(1./(N/D)**.5)*( (DN*D-0D*N)/D**2) 

I F ( X 1  .EQ.  0.0)  GO  TO  50 
SUM*SUM+(Xl+X2)/2. 


50  X1*X2 

Y1*Y1+H 
LOO  CONTI  ML' E 
IT*-H*SUM 

701*2. *PI-*RT**2*Cl*IT 

RETURN 

EMO 
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SUBROUTINE  REGEM  (T2 ,T3 ,T5 ,76 ,P ,3 ,MB  ,?R ,jPPR£G,DHDREG,7CL ! 


REAL  MB 
T-TS/2+T5/2 


p].3p-»pq 

SIGMA* ASIGMA(T.Pl) 

:r(T2  .EQ.  T3)  GO  TO  99 
DTLM*T3 -T2 

GAMMA* L . / (5  303 . *° 1 ) ** . 4 1* ( . 0  2  54*0 ) **1 . i 3/ ( . 5 6*0  TIM ) **! . A 1 
DH1*C?2( T5)-C?2(To) 


P 1-3. 14159 

AREG*P 1/4. *51 GMA*GAMMA* i MB* . A  5  4 ) * ( D  Hi *2  3  2  5 . • **1 . A 1* ( 1 . -0  r 0  RE  G ) 
***2.41*(  l.-H3R**(*2)*0HDR£G**  (-3)  )**.41*25.3 
VOL*AREG*DPPR£G**(- .41) 


RETURN 


END 


SUBROUTINE  REC'JP  ( T2H  ,T3H  ,T5H  ,TSH  ,PRH ,OPPH,D  , MA  ,MF  .ORPREC  ,VCL) 

REAL  MA.MF 

T*T5H/2.+T6H/2. 

?C*1./ (l.-OPPH) 

P=14. 7*PRH 
S I  GMA*AS  I GMA  ( T ,  ? ) 

DTLM»T3H-T2H 

GAMMA*(6303. *P )**(-. 41) *(.0254*0)**!. 13* (.56*0TLM)**( -l.il) 
QH1=C?2(T5H )-C?2(T6H) 

P  1*3.14159 
OHOREC*  .65 
CM*MA+MF 

AREC*P 1/4. *SI GMA*GAMMA* ( CM* . 4  54) * ( D  Hl*2  326 . ) **1 . 4 1* ( 1 . ^0  HO  RE  C ) 
***2 .41* (1.+0HDREC** (-3) *PC**2)** .41*35  .3 
VOL*AREC*OPPREC**(- .41) 

RETURN 

END 


SUBROUTINE  COOLER  (T7 ,T1 ,P ,D ,NTP,M3,DHDCLR,DPPCLR,VCL) 

REAL  NPP  ,NTP,MB 

IF(T7  .LT.  Tl)  GO  TO  99 

TA-T7/2.+T1/2. 

SIGMA*ASIGMA(TA,P) 

TC1*530. 

TC2*540. 

QTLM*(T1-TC:-T7+TC2)/L0G(  (T1-TC1  )/(T7-TC2) ) 

GAMMA* (6803. *?)**(-. 41)*( .0254*0)**! .18*( .56*0TIM)**  (-1  .-U ) 
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0  HI -C?  2 ( T7 ) -C?  2(11) 

PI-3.14159 

TV-TC1/2.+TC2/2. 

RK-AK(7A,P)/l.£-4 

RMU-AMU(TA)/5.4E-4 

RM*(T7 -T1)*ACP(TA)/' (1  .*10. ) 

RPR»APR(TA,?)/6.4 
NPP- .30 

RRH0-ARH0(TA,P)/64.0 
RSV=ARH0(T7 ,? )  /64.0 
ALPHA-RK*RMU** ( - . 3) *RM** .8*RPR3* .4 
CB=RM**(-2.75)*RMU**(-.25)*RRHC*RS7/NT?,'NPP 
ACLR*PI/4.*SIG»^A*6AMMA*(MB'*.454)*{0hi*232o.  ),r,r1.41*{  i.-CHCCw.RI 
**(l.+4LPHA*DHDCLR)*,r1.4i*(l.*:3'‘rQH0CLR**(-3)  )**.4i*35.3 
VCL»ACLR*0??CLR**(-.41) 

99  RETURN 
■NO 


SUB  ROUTINE  C OMBU STOR  { T3 , T4  f  T4H . T5H ,  ? R ,3 , M3 ,? R H , MA 

REAL  MB  ,MA,MF  ,M 

COMMON  /COMB/DPPRG,OPPRC ,DPPHE, OHCHR, 0PHE2 

T-T3/2+T4/2 

P»WP*(PR-OPPRG) 


,MF  ,701) 


SIGMA-ASIGMA(T,P) 

DTLM*(T5H-T3-T4H+T4)/L0G( (T5H-T3) /(T4H-T4) ) 

GAMMA- (680  3.  *?)*■*( -.4 1)*(  .0254*0)**1 .1S*(  .56*0TLM)**(-1 .41) 
DH1*CP2(  T4)  -CP2(  T3 ) 

PI-3. 14159 


TH-T4H/2+T5H/2 
TRH- ( 1 . -OPPHE ) *PRH 
PH* 14. 7* (PRH-QPPRC) 
RK-AK(T,?)/AK(TH,PH) 
RMU-AMU ( T) /AMU ( TH) 
M-MA+MF 


RM-MB/M 

RPR-APR(T,P)/APR(TH,PH) 

AL?HA*RK*RMU**(-.3)*RM**.3*RPR**.4 

AHTR-PI/4.*SIGMA*GAMMA*(M8’*' .454)*(DHl*2326,  )**1 .41*  ( 1  ,+OHOHR)  * 
*  (l.+ALPHA*0HDHR)*n.41*0HDHR**(-l.23)*35.3 
V0L*AHTR*DPHE2**(-.41) 

RETURN 

ENO 


SUBROUTINE  WHBVOL  (MB  ,MSfT6  ,T7  .MUHB  ,H2S  ,H3S  .DPPWHB  ,VOL) 

REAL  MB-, MS  ,NVJHB  ,NTU 

T-T6/2+T7/2 

CP-ACP(T) 

CPS-2.0 

CR«MS*C?S/(M8*CP) 

IF(CR  .GT.  1.)  GO  TO  10 
CMIN«MS*C?S 
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GO  TO  20 
10  CR-l./CR 
CMIM«CP*MB 
20  PS=4.Q 

,  ERS ( NWHB*CR- 1 . )  /  { MWHB- 1 . ) 

M7U*(-PS)<L0G(l.+l./CR*L0G(l.-CR*(l.-£R”il./?S))/(CR-ER^ 

*PS)))) 

0H*H3S-H2S 

VOL*  .0G045*CMlM,rMTU*MS*0H*0PP'VHB*'r(  -  .41) 

RETURN 

END 

SUBROUTINE  CGNDVOL  (MS.X.VOL) 

REAL  MS 

V0L«MS*X'*27.3 

RETURN 

ENO 
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uo 


SUBROUTINE  STMTAB 
DIMENSION  SA7LIQ ( 
OMMON  /STM/PRESS 
ATA  (COEF-TC),: 
*,  .S9573391E+03, 

* ,  -  .195136325+00 
*,- .506253495+00, 

+  ,  .555929425-03 , 
*,  .26548914E-04, 
*,  .307069195+00 , 
*,  .450CCCCCE+03, 
DATA  (CCEFpT(I),: 
* , - .323931735-14, - 
*,  .31536341E-02 , 
*,-.54272696E-2i, 

* , - .542515642-10 , 
*,  .339528975+01, 
",  .92756543E-17,- 
",  .399265985-08,- 
OATA  (CCSFFTC I )  ,1 
*,-.lA722587E-21, 
",  .274555055-12,- 
*,  .985127545*00, 
*,  .72857 180E-22 , 
*,-.70679233E-14, 
\  .30734203E-Q4, 
*,-.335080325-23, 
DATA  tCOErrTII)  ,1 
*,  .1558777QE-07 , 
",  .550000005*03, 
*,  .449796725-12, - 
*,  .12399614E-02, 
\  .129453715-21, - 
*,  .17548Q01E-12, 
*,  .101509745+01, 
DATA  (COEFFT(I).I 
*,  .592173555-16,- 
*,  .595984105-03, 
\  .250000005+03, 
*,  .49239634E-17, 
*,-.3£546496E-06, 
*,  .190000005+01,- 
*,-.187961095+03, 
OATA  (COEFFTdM 
*,  .111994905+04, 
*,  .62379937E+05 , 
*,  .35122466E-04,- 
*,  .30000COQE+0  1 , 
*,  .52122736E+03,- 
*,  .12340954E-04, 
*,  .18000000E+01, 


( N I ,  X 1 , 7 1 ) 

3) .STEAM (3) .COEFFT  !227  ) 

,te:-p  ,  dual,  satliq,  steam 

-1 ,30)  /  .2208222I5+Q4, 

.6046426 1E+01,  .7CCC00GCE+Q1, 


.522799525-00 
.725957205+00 
.2000000  OE’-Ol 
.91349306E-04 
.262413435+01 
.5SC00C00E+C3 
=  31,  60) 

.1 1734466E-11 
.69510255  E— 0 1 
. 1 27042665 -IS 
.447973225-07 
.A3Q1Q1S3E-02 
.954706565-15 
.99567212E-C6 
=  51,  90) 
.56553573E-19 
.133353025-10 
.25CC0CCCE+02 
.332255065-19 
.901142015-11 
.2175269 15-0 1 
.353197925-20 
=  91*120)'  "  . 
.3126444GE-05 
.947431475-13 
.51337380E-10 
.1279855 15^0  L 
.553731295-19 
.23202424E-09 
.21848036E+03 
=121,150) 
.183453225-14 
. 4557591 5E-06 
.50084539E-23 
.15250351E-13 
.142857095-02 
.110850775*05 
.156997975+03 
=151,130) 
.150000005*01 
.29474498E*05 
.1234Q832E+04 
.130000005-04 
• . 1 1 6829475 -0 1 
.135551265-03 
.2282209 1E+02 


,  .311489535-01 
,  .1 458569-5-02 
.-.7413344CE-06 
,  .156799225-02 
,  .4496061 5E-K)2 
,-.987368385-22 
/  .200411395-17 
,  .H244476E-07 
,  .250000005+02 
,-.569751635-16 
,-.231047555-04 
,  .549273775-02 
,-.290046775-12 
,-.591579925-03 
/  .550000005*03 
,  .477256055-17 
,- .311603335-06 
,  .470000005+03 
,  .300755755-18 
,  .581496425-09 
,  .250000005+03 
,  .25729265E-13 
/  .5C5C6309E-13 
,  .170025395-01 
,  .24971473E-15 
,  .491789865-08 
,  .549273775*03 
1381757 4E-17 
,  .361459675-06 
,  .45000000E-02 
/-  .211344385-20 
,-.66417453E-12 
,  .12325099E-02 
,-.156590705-21 
.7  40683205- L  2 
:,  .36 7534Q5E*00 
,-.159945875*05 
,-.202133425*03 
/  .150404555*03 
,-.148635855*06 
,-.139075365+05 
',  -.539298275 ■'■02 
,  .708489905*01 
, -.12001338E-00 
,-.150993195-02 
.187766495*03, 


.135526805-00 

.513426325-01 

.124042195-03 

.732248225-06 

.250325205-01 

.179267225-03 

.192012305-19/ 

-.421126595-15 

-.'24196525-05 

,2C5"55485-22 


-.1244A7’'?5-i9 
.252522925-10 
.731422745-00/ 
.21343036E*03 
-.533491625-14 
-.119563545-03 
.550CCC0CE-03 
-.728986515-15 
.367129815-07 
-.734591295-25 
-.315919045-15/ 
.136845995-10 
.45CGCCCC5-C3 
-.383025385-14 
.483253955-05 
.25C0CQQ0E*03 
.625201065-14 
.125064615-03 
.550000005*03/ 
.1532703SE-13 
.145895635-10 
.749697975+00 
-.460991195-13 
•10479899E-08 
.1 'COCOOOE+Ol 
.914076885+04 
-.486501755*02/ 
-  .130505475+03 
-.13414819E*06 
-.5785387  55-02 
.120268265+04 
.16507788E-02 
.204119915-01 
.20CQC0005-01 
.320805045-02/ 
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DATA  (CCEFFT(I),  1*131, 210)  /- .3ZU14G8E+01, 

*,  . 1540083 58+04,  .137466578-02,  .5397387 38+02, 
*,  .13000CQ0E+Q4,  .130C0Q0QE+01,-.542698aCE-05, 
*,  .28879248E+01,  .11Q34G6GE-05 ,  .472250018-03, 
*,  .69416973E-05, -.50892120E-02,  .1G9S19ZZE+01, 
*,  .13000C0QE+G1 ,- .331Q1208E-G3 ,  .215519158+01, 
*,-.77236ol8£-Q3,  .3675365CE+00, - .2545012IE+C3, 
*,-.164271068+01,  ..S1477766E+03,  0.0, 

DATA  (COEFFT(I), 1=211,237)  /  0.0,  G.O, 

*,  Q. 0,0. 0,0. 0,0 .0,0. 0,0.0 ,0.0,0. 0,0. 0,0. 0,0.0, 
*,  . 1 4000000E+0 1 , - . 1 5 47 9259E -04,  .193S1S56E-01, 
*,  .23329644E-04,  .114937918-02,-. 131392088*02, 
*,-.125923998-01,  .5510115CE+01/ 

M-NI 
X-XI 
y=y  i 
I-N-3 
LINK-0 
PRESS* X 
DELTA1=0 .0 
ERR0R1-0.0 
8RR0R2-0 .0 

IF  (M-o)  1000,1000,1030 
1000  IF  (I)  1020,1020,1010 
1010  LINK-1 
N-I 

1020  IF  (X-COEFFT(N) )  1040,1040,1520 
1030  N-I 


.424491018-0 

.207198418-0 

.964138028-0 

-.941136008-0 

.13CCCCQCE-0 

.533002918-0 

.252497078-0 

0.0 

'4*0)  ^  j  '■J  • 

.9C000CG08-0 

.245939098-0 

.588602138-0 


1040  GO  TO 
1050  X*Z 


(1 130, 1G60, 1190, 1230, 1080,1030, 1180, 1330)  ,.‘j 


N=2 

1060  TEMP-X 

GO  TO  1190 
1070  PRESS=EXP(Z) 
GO  TO  1150 
1080  11=2 
JJ-3 

TOLER*  .005 
GQ  TO  1100 
1090  1 1-3 


JJ=2 

TOLER- .000005 
1100  H-Y 
1110  LINK-2 
1120  N-I 

1130  X=ALOG( PRESS) 
PLOOX 
GO  TO  1190 
1140  T8^=Z 
1150  X-T8MP 
>1 
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M*4 

30  TO  1190 

UoO  satl:q(j)*z 
J*0+1 
M*N+1 

IF  (0-3)  1190,1130,1170 
1170  IF  (LINK)  1640,15-10,1300 
1180  H*X 
S*Y 

x*s 

1190  i=23*M-19 
K»1 

IF  (X-COEFrT(I))  1200,1200,1210 
1200  .<-12 
1210  I=I-< 

ARS»X-CaeFFT{ I) 

7-Q^Q 

00  1220  <*1,10 

1*1  +  1 

1220  Z*Z*ARG+CQEFFT( I ) 

30  TO  (1140,1070,1050,1150,1160,1150,1410) ,N 
1230  TEMP=Y 

30  TO  1300 

1240  £RROR*H -STEAM ( 1 1 ) 

IF  (A8S( ERROR) -TOLER)  1540,1250,1250 
1250  IF  (ERRQR1)  1260,1280,1260 
1250  £RR0R1*ERRQR1 -ERROR 

IF  (ASS (ERRCR1) -TOLER)  1290,1290,1270 
127Q  SLOPED ELTA1 /ERR0R1 
1230  0EITA1*3L0PE*ERRCR 
ERR0R1*£RR0R 
1290  TEMP  *TEMP  +0 EL T  A1 
1300  IF  (PRESS-COEfFT(l)i  1320,1320,1310 

1310  IF  (( .0Q00n78*9RESS+.39411)-,'PRESS+382. 1-TEMP)  1320,1320,1529 
1320  T=*  .55555556*TEf-P+255. 38223 
TLQG*ALOG(T) 

TAIM.O/T 

COEFFT  (222 )  *TAU*TAU*136210 .06 
C0EFFT(217)«EXP(C0EFFT(222)+7.8791476-TL0G) 

COEFFT (221) *162460. *TAU 
C0EFFT(220)»125970.*TAU 
C0EFFT{2I9)*£XP(149. 27765-23. *TIQG) 
C0EFFT(212)*1.39-C0EFF7(217) 

COEFFT (21 1 )*82 .546 -COEFFT (22 1 ) 
C0EFFT(21Q)*.2182S*T-CQEFFT(220) 

COEFFT  (209  KOEFFT  (219 )-  .000363  5*T 
C0EFFT(222)*(C0EFF7(222)+C0£FFT(222)+l.Q) *C0EFFT(2 17) 

COEFFT (21 7 ) *COEFFT (212 ) -COEFFT (222 ) 

C0EFFT(214) *C0EFFT(217) /COEFFT (2 12) 

COEFFT  (216)  *CQEFFT  (211)  COEFFT  (21 4 )+4 1 .272  -COEFFT  (22 1 ) 
CQEFFT(215)*C0EFFT(2105*C0EFFT(214)-.5*C0EFFT(22Q5 
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C0EFFT(214)3C0EFFT(209)*C0EFFT(214)+1.8461538*C0EFFT(219) 

COEFFT (22 1 ) 3 .50+C0EFFT (211) -COEFFT (21 5 ) 

COEFFT(22Q)3 .25*CG£FFT(210)-CQEFFT(2 15} 

CCEFFT (21 9 ) 3 ,O75923O77*C0EFFT (209 ) -COEFFT (214) 

C0EFFT[213)3(  .0037299965*T+14. 531313) *T+17906. 513+472. 2493*rLCG 
COEFFT (223 )3- . 007459993 1*T-14. 53131 3*TL0G+ 2 . 4524439*47 2. 2493 7 
**TAU 

1330  P* .Q53046190*PRESS 
PLCG*ALOG(  PRESS) 

C0EFFT(213)3C0EFFT (212)*TAU 
C0£FFT(226)=C3EFFT(213) 


COEFFT (225 ) 3CC£FFT (226 ) *C0£FFT (225 ) 
COErFT(224)3COEFFT(225)*COEFFT(225) 

COEFFT  (224X0EFFT  (224)*C0EFFT  (224)  -COEFFT  (226  ) 
C0£FFT(225) =CC£FFT(225) *CCEFFT(213) 

COEFFT  (213)SP*TAU 


STEAM( 1 )3( ( ( COEFFT { 209 ) ’COEFFT ( 224 )+COEFFT( 2 10) }*CCEF 
*CO£FFT (211 ) )*COEFFT (225  )n)0EFFT (21 2 )+4 .55504, 'COEFFT (21 
**.0150135 

STEAM{2)*(  ( ( (C0EFFT(214)*CC£FFT(224)+C0£FFT(215)  )*COE.; 
*CQEFFT (216) )*COEFFT (225 )+CGEFFT (217)) *P*CO£FFT (218 ) ) * . 
STEAM(3 )3( ( ( (C0EFFT(219)*CGEFFT(224)+CC£FFT(220) )*CCEF 
*C0EFFT (221) )*COEFFT (225/+CCEFFT (222) !*CC£FFT (213 )+CO£F 
*+4.55504*PLGG)*(-. 02419825) 

QUAL31.0 

I=LINK+1 

GO  TO  (1640,1640,1340,1240,1340,1450),: 

1340  X*STEAM( ! 1 ) 

IF  (H-X)  1350,1640,1360 
1350  Y-SATLig(II) 

qUAL=(H-Y  )  /  (X-Y  ) 

X*1 .O-QUAL 

STEAM(1 )3CUAL*STEAM(1)+X*SATLIQ(1 ) 

STEAM (0  J ) -Q U AL *STEAM ( J J ) +X*S ATI I Q ( J  J ) 

STEAM(II)3H 

IF  (LIMK-2)  1640,1640,1460 
1360  I3( I I-l)*ll+154 
J*1 


-'(225) 


225)- 

35563 

(225)- 

(223) 


X*PLOG 

1370  00  1380  K*2Q9 ,219 

COEFFT(K)3COEFFT(I) 

1380  1*1+1 

X*X-C0EFFT{209) 

Y3Y-C0EFFT(210) 

CQEFFT(2Q9)*(CQEFFT{21i)*X+CQEFFT(212) )*X+C0EFFT(2 13) 
COEFFT (210 ) -( COEFFT (21 4) +X+COEFFT (215)) *X+C0EFFT (216) 
C0EFFT(211)3( (COEFFT(217)*X+COEFFT(218) )*X+CQEFFT(219) )*Y 
Z*COEFFT (211 )+COEFFT (210) 

SLOPE*C  OEFFT ( 2 1 1 ) +Z 
Z-Z*Y+C0EFFT(209) 

GO  TO  (1390,1430,1530,1540) ,J 
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1390  LIMK=3 

IF  (Z-TE'-’P)  1200,1300,1400 
1400  TSMP=Z 

GO  TO  1300 
1410  J-2 
1=227 
:<=h 

IF  (Z-H)  1420,1370,1370 
1420  1=137 

GO  TO  1370 
1430  DPOS=SLOP£ 

?RSSS=£XP{ Z) 

PLOG*Z 

IF  (1-227)  1520,1520,1440 
1440  EIMK*4 
11=2 
30=3 

GO  TO  1120 
1450  DH=H-ST£AM(2) 

STE  AM(  3 )  =ST£A.‘U  3 )  *0  H  /  ( T  Zl-V +459.69) 

1460  QS=S-ST£AM(3) 

DELTA1=STEAM(3 )-£RRORl 
IF  (ERR0R1)  1470,1490,1470 
1470  IF  (A8S(Q£LTA1 )- .0QG005)  1500,1500,1480 
1480  DPDS=(PLQG-PL0G1) /DELTA! 

1490  PLGG1=PL0G 

ERR0R1=STEAM(3) 

1500  PRESS=PRESS*(1 .Q+OS*OPOS) 

IF  (LIMK-4)  1510,1510,1550 
1510  IF  ( ABS(QS)- ,QCGCQ5)  1(40,1120,1120 
1520  LINK=5 
X=H 
Y=S 
1=198 
J=3 

GO  TO  1370 
1530  TEMP=Z 

DTDS=SLOPE 

X=PLOG 

Y=H 

1=165 

J=4 

GO  TO  1370 
1540  OTDH-SLOPE 
GO  TO  1300 

1550  TEM>=TEM>+OS*OTDS 

IF  (A8S(DH)-.005)  1610,1560,1560 
1560  3£LTA2*Q  H-ERR0R2 

IF  (ERR0R2)  1570,1590,1570 
1570  IF  {A8SOELTA2 )-  .005)  1600,1600,1580 
158Q  OTO  H*  ( Ttf'P  1  -TEMP )  /DEL  T A2 
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1590  TEMP1«TEMP 
ERRCR2=0H 

1600  TEMP*TEMP+QH*9TDH 
GO  TO  1320 

1610  IF  (A8S  (OS )-  .000005)  16*10,1500,1500 
1520  WRITE  (6,1530)  MI.XI.YI 
1630  FORMAT  ( 13 .2E14.3 ,13'H  OUT  OF  RANGE) 
STOP 

16*10  RETURN 
END 
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L 


2.00 


12  T2  '  V*^nM  T  . 

n  ^  i  *  w 


P.AT-C  s 


WORKING  PRESS  "JRE  =  1  35 . 0 
OCTAL  PRESS  USE  DRC ?=  . : C 
MASS  FLOW  (L cm/s)  =*±44. 3? 


STATE 

TEMP (deg  R) 

5  Ur.DlAR  Y :  THER  MAL  E?: 

< 

** )  •  ,■** 

.  n 

HEATER  EPF 

2 

675.7 

3FC ( Lbm/K? - 

3 

1663.3 

7pr  ’  )' 

*  sj  x.  V 

*r 

2C60 . 0 

POWER (H?) =2 

-/ 

13C4.2 

WORK  FRAC  2' 

6 

311.2 

?? 

734. 0 

7 

779.3 

COMPONENT 

IP?  I HD 

VCL( cu  ft) 

REOEN 

• C7S  • 7^^* 

15-5  EFF=0 . 33 

TURBINE 

1 0  .  C  EFF  ( ? ) =C . -1 

wwC  w-T.  rX 

.00?  .250 

159.3 

^  t-^o  rcc^o 

VSa  LH  A  AwO 

29.3  EFF;?) =0.33 

HEATER 

—  r  r  z 

•  w  0  \j  y 

3C ILER 

.014 

T  IT  ^  O  ^ « 1 

— »  r  -  C  . 

ITT  \  m^rp  VP-  T  . 

.  — A--A  • 

PRESSURE  RA 

TTC  =  -.79 

TOTAL  PRESS UR 

E  DROP* .43 

AIR  FLOW  (Lira/ 

'3)=c7.34 

FUEL  FLOW (Ida 

/3)=  3.?7 

STATS 

TEMP( deg  R) 

Ih 

535.0 

2h 

335.2 

1727.- 

5054.5 

5h 

1763.3 

6h 

945. 0 

7h 

760.  C 

COMPONENT 

DP?  DKD 

7CL( cu  ft) 

RECUP 

.321  .65c 

2.5  E?F  =  G  .96 

TURBINE 

1.1 

HEATER 

.012  .310 

67.4 

COMPRESSOR 

16.0  EFF(?)=0.35 

RANKINS  CYCLE 

i  SATURATION 

PRESSURE*  56.7 

STEAM  FLOW ( L5m/ s ) =  2.96 

FLOW  RATIO =0. 

0067 

STATE 

TSM?(deg  R) 

PRESSURE (csi) 

Is 

540.0 

0.5 

2s 

555.0 

56.7 

SF 

749.0 

56.7 

3s 

776.2 

56.7 

4s 

540.0 

0.5  Q UAL =0.544 

COMPONENT 

VCL(cu  ft) 

30 ILER 

162.3 

TUR3INE 

4.3  E?F( IS ) =C . 38 

CONDENSER 

63 .3 

^  *  ?  0  V 


r  7  r>  n  -> 

1  ^T»u  w 


30 


i  /\  v/i 


mvA  **w*i  ^  -  -  —  • 

PRESSURE 

WORKING  PRESSURE  =135-0 

TOTAL  PRESSURE 

DROP*. 10 

.'lAsS  FLC.i  (~cr./ 

s)*2S7.63 

T^  \  -T» 

EM?( deg  R) 

l 

340 . 0 

2 

769.3 

3 

1538.3 

4 

2060. 0 

5 

1643 . 6 

6 

374.7 

?? 

341.2 

7 

3  31 .3 

COMPONENT 

DP?  3 HD 

REGEN 

<\  £  r\  2  7' 

.  -  >  ^  ~  J 

turbine: 

TC 

w  v.  V  l—*  - 

.024  .35 

COMPRESSOR 

HEATER 

.0004 

BOILER 

.026 

HEATER  CYCLE; 

PRESSURE  R 

RA--C 


-  *J 


SUMMARY :  THERMAL  E? ?  *0  .  4 

V,-^"T3  TTT  s  Q  t  C 

3  F  0  ( L  cn/H?  -  r*r )  = 
7CL  TOTAL (  cu  ft 
POWER  C  H? )  =  2  5  0  C  C 
,jCr:<  frac  by  sr: 


C^-9C 


*  *  t  _  ,  ,  «i 

^ ^  v  C  U  *  - 


-T*.*  > 

TTT?-n 


50 
38 . 9 
20.5 


-- (  -r>  \  -•  n 

ir -  \ r  ;  -  v 


0 . 3; 

EF?=C .  ?4 


»  ,fl  —r<  -  !<  7Q 

r\-  «V»  “■  •  (  ; 


:rc?=.^3 


TOTAL  PRESSURE  - 

ALR  FLCW'(Lcn/s}  =  40.99 
FUEL  FLOW  ( L’am/s )  =  3-^G 

STATE  TEMP ( deg  R) 


Ih 

2h 

3ii 

Lh 

5R 

6h 

7h 

COMPONENT 
RE  CUP 
TURBINE 
HEATER 
COMPRESSOR 


535-0 
335.2 
1336 . 5 

4932.9 

I63S. 3 
9*5.0 


RANKINE  CYCLE; 

STEAM  FLOW ( Ltm/s ) =  2.43 
FLOW  RATI0=0.0086 


DPP 

DKD 

7CL( cu  ft) 

.310 

.650 

2.2  EFF=0 .93 
0.9  EFF ( ?) =0 

.022 

•  358 

23 . 6 

12*7  EFF( P) =0 

SATURATION 

PRESSURE =  128.0 

STATE  TEMP (deg  R) 

Is  540.0 

2s  555.0 

SP  306.2 

3s  339-7 

4s  540.0 

COMPONENT  VOL( cu  ft) 

BOILER  38.3 

TURBINE  3-*  SFF(  IS )  =0 . 38 

CONDENSER  55 *5 


PRESSURE (psi) 

0.5 

128.0 

123.0 

128.0 

0.5  QUAL=0 . 319 
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BRAYTCN  CYCLE  i  PRESSURE  RATIO*  -.  CO 
.<CRKINu  PRESSURE  *  135.0 
TOTAL  PRESSURE  DROP*.  10 
MASS  PLOW  (Lbm/3)=2  39.?1 


cm  '  me  me 
O  ~  ~  — 

01? ( deg  R) 

SUMMARY :  THERMAL  E??*C  .  4 -a  7 

l 

5-C .  0 

HEATER  EPP  =  0 . 9  3 3  ( PHI  =  C . 9  0 ) 

2 

343.9 

3PC  ( Lbns/H?  -hr )  =  C  .  - 7  3 

3 

14=3.3 

70 L  TC^-L-'  :u  '5£2.2 

4 

ZOoQ.O 

power! h?7 =25000.0 

5 

153^.9 

WORK  PP.AC  BY  STM  CYCLIC.  Cc  3 

/ 

0 

927.0 

876.3 

( 

860.6 

f+  ^  ''5,^  rTC  Y  ’ m 
u  w.^ii  — 

-~-_o 

VCL( cu  ft: 

REuEM 

.021  .93C 

13.2  EPP  *0 . 3c 

T  UR  3 1  ME 

^  T"?7  l  \  q  * 

-  r  rs 

w  w  v*  L-  •  -*  IY 

.032  .44A 

t?  •  t 

CCMPRESSCR 

13.4  EP?( ? ) =0 . 33 

HEATER 

.0003 

BOILER 

.0  46 

EP?=C  .9^ 

HEATER  CYCLE : 

PRESSURE  RAT 

10*  4.^9 

TOTAL  PRESSURE 

DROP*. 43 

AIR  PLOW (L bn/s; 

=39.59 

PUEL  PLOW ( Lbm/s ) =  3.29 

stats  te 

IMP  (deg  R) 

1  h 

535-0 

Zh 

a  -jc  0 
•  - 

3h 

1494.5 

Ah 

4355  •  0 

5b 

1553. 3 

6h 

945.O 

7h 

760.0 

COMPONENT 

DPP  DHD 

7CL(cu  ft) 

RECU? 

.299  .650 

2.1  ZP?  *0.92 

TUR3IME 

0.3  Z??(?)=0.39 

HEATER 

.033  .406 

17.7 

COMPRESSOR 

12.0  E?P( ?) =0 .35 

RANKINS  CYCLE: 

SATURATION 

PRESSURE*  199.9 

S TEAM  PLOW ( L  cm/s )  =  3-19 

PLOW  RAT 10*0.0133 


STATE 

TEMP ( deg  R) 

PRESSURE (?si) 

Is 

540.0 

0-5 

2s 

555.0 

199.9 

SF 

341.3 

199.9 

3s 

392.0 

199.9 

4s 

540.0 

0.5  0UAL=C . 310 

COMPONENT 

VOL(cu  ft) 

30  -i.ER 

119.0 

TUR3INE 

5-1  EPP 

‘(IS)  *0.33 

CONDENSER 

70.5 
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3RAYTCN  CYCLE :  ??.£ 
WORKING  PRESSURE =  13 
TOTAL  PRESSURE  DROP 
MASS  FLOW  (L cm/s)  =  2 


COMPONENT 
REOEN 
T  UR  3  INS 
COOLER 
COMPRESSOR 
HEATER 
3CILER 


2060.0 

1-58.1 

972.1 

3?6.8 

3L6.9 

-  2; 


r^lCU  I 

.  < 

*2 !  7 
107.3 
16.7 


.  “HE?  UAL 
HEATER  EEP=0. 
S?C(Lbm/K?-hr 
7CL  TOTAL  Uu 
POWER  (HP  ’=25C- 
WORK  “RAC  BY  . 


E??=0 .33 


=0 .-cC 

933  <E: 

) =C .-6 


5CCC.C 


0003 

068 


E7?(  ? )  =( 


E??=0 . 9^ 


HEATER  CYCLE j  PRESSURE  RATIO*  4 . 79 
TOTAL  PRESSURE  CROP*. 43 
AIR  PLOW ( L'cm/s )  =  33.-9 
FUEL  FLOW ( Lbm/s ) =  3.19 


'ATE 

TEMP (deg  R) 

lh 

535-0 

2h 

335-2 

3b 

1^27.- 

4h 

-798.9 

5h 

1491.3 

6h 

945.0 

7b 

760.0 

COMPONENT 

-'■D'O 

-fc.  A 

nur) 

wr. u 

7CL( cu  f 

*  / 

RECU? 

.239 

.650 

2.1 

EF?*C .90 

TURBINE 

0.3 

EFF(  ?)  =0 . 89 

HEATER 

.  044 

.454 

12  .4 

COMPRESSOR 

U.5 

SF?( ?) *C . 35 

RANKLE:  CYCLE:  SAT  UR  AT 
STEAM  FLOW (L'om/s)*  5.21 
FLOW  RATI0=0.02L5 


SATURATION  PRESSURE*  200.0 


STATE  TEMP (deg  R)  PRESSURE ( osi) 

Is  540.0  0.5 

2s  555.0  200.0 

SF  341.3  200.0 

3s  937.1  200.0 

4s  54o.o  0.5  QU 

COMPONENT  7CL(cu  ft) 

3CILER  298.1 

TURBINE  10,9  SFF( IS ) =0 . 88 

CCNEENSER  117.4 


QUAL*0 . 826 


0\0 


r 


3RAYTCN  CYCLE:  PRESSURE  RATIO  =6.00 
WORKING  PRESSURE*! 35 . 0 
TOTAL  PRESSURE  DRC?=.tO 
MASS  FLOW ( Lem/ s )  =  1 9 7 . 0 0 


state 

TEMP (dag  R) 

SUMMARY :  THERMAL  EFF* 

1 

<40.0 

HEATER  EFF=0 

2 

959> 

3  F  C  ( L  'cm/H?  -  hr 

3 

1343.3 

YCL  TOTALicu 

/« 

2060.0 

POWER (HP ) =25C 

5 

1396.2 

WORK  FRAC  BY  : 

6 

1011.3 

?? 

376. 3 

7 

335.3 

COMPONENT 

DP?  DKD 

7CL( cu  ft) 

REGEN 

.010  .950 

14.4  EF?=0 . 38 

TURBINE 

2 . 2  EF? ( ? ) =  C  .  9 1 

GCGLER 

.016  .623 

129.7 

COMPRESSOR 

15.3  E??(?}=C.3c 

HEATER 

.0002 

BOILER 

.074 

EF?*C  .94 

HEATER  CYCLE: 

PRESSURE  RATIO*  4.79 

TOTAL  PRESSURE  DROP*. 43 

AIR  FLOW  ( L’om/s )  =  3 S .  i  7 

FUEL  FLOWCL-o/ s ) -  3-17 

STATE 

TEMP( deg  R) 

lh 

535-0 

2h 

335.2 

3h 

1375.4 

4h 

4755.3 

5h 

1443.3 

6h 

945.0 

7h 

760.0 

COMPONENT 

DPP  DHD 

VOL(cu  ft) 

RSCU? 

.273  .650 

2.1  EFF=0 . 69 

TURBINE 

0.8  E?F(?)=0.39 

HEATER 

.055  .502 

9-5 

COMPRESSOR 

11.4  EFF( ?) =0 . 35 

RANKINE  CYCLE:  SATURATION 

PRESSURE*  200.0 

STEAM  FL0W(L' 

bm/s)=  6.70 

FLOW  RATIO =0 

.0340 

STATE 

TEMP( deg  R) 

PRESSURE ( os i) 

Is 

54  0,0 

0.5 

2s 

555.0 

200.0 

S? 

341.3 

200.0 

3s 

976.3 

200.0 

4s 

540.0 

0.5  QUAL=Q . 333 

COMPONENT 

V0L( cu  ft) 

30ILSR 

523.6 

TURBINE 

16.2  EFF( IS  5  =0 . 38 

CONDENSER 

1530 

34 


i 


l/l  C»  1 1. 


3RAYTCN  CYCLE;  PRESSURE  RATIO*?. CO 
WORKING  PRESSURE =135 . 0 
TOTAL  PRESSURE  DRC?=.10 
MASS  ?low(L'em/s)  =137  *43 


CiT>  * 

?ESS?(  deg  R) 

SUMMARY s THERMAL  E??  =  0 .46  3 

1 

540 .0 

HEATER  E??= 0.933  ( 

2 

10C7.0 

S  FC ( Lbm/H? =hr ) * C . 4  5 

3 

1304.9 

70 L  TOTAL ( cu  ft) =11 

4 

2G60.0 

POWER ( HP) =25000.0 

5 

l>5-5 

WORK  FRAC  BY  STM  CY 

6 

1047.6 

pt> 

376.3 

7* 

325.2 

COMPONENT 

DP?  SHE 

7CL(cu  ft) 

REGEN 

.CIO  .950 

13.6  I??=C . 33 

TURBINE 

1.9  2??(?)»0.91 

COOLER 

n  -  -5  70  /*> 

•  f  -  w 

155.9 

COMPRESSOR 

150  EF?(  ?)  =0 . 33 

HEATER 

.0CC4 

30  t  j.R 

.073 

S??=0.9 4 

•-nr'rf'tro  fv^rr. 

PRESSURE  RA 

TIC  =  4.79 

TOTAL  PRESSOR 

E  IRC ?= .-3 

AIR  ?L0W(Ltm/s)=  33-22 

FUEL  FLOW { Lom/s ) =  3-17 

STATE 

TEMP (deg  a) 

Ih 

535.0 

2h 

835.2 

3h 

1333.3 

4h 

4721 .2 

5h 

1404.9 

oh 

945.0 

7h 

760.0 

COMPONENT 

3P?  DKD 

VCL(cu  ft) 

RECUP 

,267  .650 

2.1  EFF=0 . 37 

TURBINE 

0.3  EF?(?)=C.39 

HEATER 

.066  .550 

7.7 

COMPRESSOR 

11.4  EFF(?) =0.35 

RANKINS  CYCLE 

SATURATION 

?RESSURE=  2C0.0 

STEAM  FLOW ( Lbm/s ) =  7-92 

FLOW  RAT 10=0. 

0422 

STATE 

TSMP(deg  R) 

PRESSURE ( os i) 

Is 

54o.o 

0.5 

2s 

555.0 

200.0 

SF 

341.8 

200.0 

3s 

1012.6 

200.0 

4s 

540.0 

0.5  QUAL=0 . 349 

COMPONENT  VOL(cu  ft) 
3CILER  782.4 


TURBINE  21.0  SF?(IS)=0.88 

CONDENSER  133.6 


’ll 


3RAYTC  N  CYCLE: 

PRESSURE  3 A 

.T 10=  3.0C 

WORKING  PRESSURE*  135-0 

TOTAL  PRESSURE 

CROP*  .10 

MASS  FLOVHLtn/s)  =  131.22 

STATE  TEM?( deg  R) 

SUMMARY  : THERMAL  E?F=C.~cC 

4 

J. 

540 . 0 

HEATER  EPF =0.9 33  v? 

2 

1049  .3 

SFC ( Lbm/K? -Hr ) =0 . 4  5 

3 

1272.5 

VP  r  m  *  7  r  r  •  \  —  *  r 

<  4^  iw.  rv-j  -  w  y  > 

4 

2060.0 

POWER ( HP ) =25000. C 

5 

1302.9 

WORK  PR AC  3Y  STM  1Y 

6 

1030.2 

PP 

376.3 

7 

316.3 

COMPONENT 

DPP  D HD 

YCL(cu  ft) 

REGEN 

.CIO  .950 

l"5.!  EPF  =  0.3  3 

TURBINE 

i!S  i?p(?)=c.9i 

COOLER 

.010  .312 

131.1 

COMPRESSOR 

15.1  EP?(?) =0.35 

HEATER 

.0004 

BOILER 

.030 

EP? =0.94 

HEATER  CYCLE: 

PRESSURE  RAT 

10 »  4.79 

TOTAL  PRESSURE 

DROP*. 43 

AIR  PLOW ( Lbm/s 

)  =33.45 

FUEL  PLOW ( Lbm/s )*  3.19 

STATE  T: 

SMP( deg  R) 

lh 

535-0 

2h 

335-2 

3h 

1293.2 

4h 

4692.5 

5h 

1372.5 

6h 

945.0 

7h 

760.0 

COMPONENT 

DPP  DKD 

VOL(cu  ft) 

RECUP 

.256  .650 

2.1  E??=0 . 36 

TURBINE 

0.3  E?P(?) =0 . 39 

HEATER 

.076  .598 

6.5 

COMPRESSOR 

11.5  EPP(?)=0 .35 

RANKINE  CYCLE: 

SATURATION 

PRESSURE *200.0 

STEAM  PLOW ( L'om/ s ) *  3.97 
FLOW  RATIO *0.0495 

STATE  T: 

SMP( deg  R) 

PRESSURE ( os i) 

Is 

540.0 

0.5* 

23 

555-0 

200.0 

SF 

341.3 

200.0 

1045.2 

200.0 

4s 

540.0 

VCL( cu  ft) 

0.5  0UAL=0 . 359 

COMPONENT 

BOILER 

1033.1 

TUR3INE 

25.5  EFF 

( IS) =0.38 

CONDENSER 

210.2 

36 


( 

l 


I 


i 


^nU.  *U 


1 


3RAYTCN  CYCLE:  PRESSURE  RATIO*  9-00 
W ORKING  PRESSURE =135-0 
TOTAL  PRESSURE  DROP*. 13 
MAS S  FLOW  ( L bm/s )  *  1 76 . 9 5 


LIP  (deg  R) 

SUMMARY: 

THERMAL  EPF=C 

540.0 

HEATER  Z??=0 

IO89.O 

SFC(Lbm/K?-ar 

1245.0 

7CL  TOTAL (cu 

2C60 . 0 

POWER (H?) =250 

1266.2 

WORK  PP.AC  BY  ; 

1110.2 

376.8 

303.2 

DP?  UHL 

7CL(cu  ft) 

.010  .950 

12.7 

SFF=C.cc 

1 .4 

S?P( ?) =c . 91 

.010  .904 

197.5 

EPF ( ? ) =C . 23 

15.0 

.  0004 

.080 

E??=0 .94 

PRESSURE  RAT 

10=  4.79 

DROP  =  .4  3 

"i'~ 

2 

-j 
-/ 

4 
< 

/* 

0 

7? 

7* 

COMPONENT 
REGEN 
TURBINE 
COOLER 
COMPRESSOR 
HEATER 
BOILER 

Tvr,,f  t 

TOTAL*  PRESS Ui 
AIR  FLCW( Lbm/s )=  33.30 
PUEL  FLOW ( Lbm/s ) =  3-22 

STATE  TEMP (deg  R) 

lh  535-0 

2h  335-2 

3h  1263.3 

Lh  4663 .2 

5h  13^5-0 

6h  9^5*0 

7h 

COMPONENT 
RECU? 

TURBINE 
hnATER 
COMPRESSOR 

RANKINS  CYCLE:  SATURATION  PRESSURE®  200. CO 
STEAM  PLOW ( Lbm/s )*  9.90 
PLOW  RATIO =0.0560 

STATS  TEM?( deg  R) 

Is  5^0.0 

2s  555-0 

S?  841.3 

3s  1075.2 

4s  540.0 

COMPONENT  7CL(su  ft) 

30ILER  1439.0 

TURBINE  29.3  EPF (IS) 

CONDENSER  234.5 


.56 

>33  v 
=  0.46' 


:0 . 90) 


STM 


0.260 


760. 

0 

t) 

DPP 

DHL 

7CL(cu  f 

.245 

.650 

2.2 

E??=0 . 35 

0.3 

EPF( ?) =0 . 89 

.087 

.646 

5.7 

EFF(?) =0.35 

11.7 

PRESSURE (?si) 
0.5 

2C0.0 
200.0 
200.0 
0.5 


!0 . 38 


£UAL=0 . 367 


/ 


37 


BRAYTCN  CYCLE:  PRESSURE  RATIC-lU.uu 
FORKING  PRESSURE =135  •  0 
rCTAL  PRESSURE  DRC?=.iO 
MASS  FLOW  ( Lem/ 3 ) =1 ?4 • CO 

3T4TS  TEMP  (deg  R)  SUMMARY: 

1  5-o.o 

2  1125.1 

t  1  ■’21  1 

2  2060 !  0. 

5  123^.2 

6  1138.1 

??  376.7 

7  800.9  . 

COMPONENT  DP?  DhD  /OiAcu 

REGEN  .CIO  .950  -2- ^ 


3 UKMAR Y : THERMAL  ZFF  2  C . 

HEATER  EFF=0.9 
S  ?C ( L bm/H? -hr )  = 
7CL  TOTAL! cu  ft 
PC'.vER  ( HP)  =2  5  COC 
.VC RX  FRAC  BY  ST 


TUR3INE 
CoC  u^.R 
COMPRESSOR 
HEATER 
BOILER 


.996 


YOL(cu  f! 

12.  4 

1-3 

2j2.3 

14.9 


.0004 

.082 


EFF=0 .38 
ZFF ( ?) =C .91 

EF? (?) =0 . 83 

EF?=0 .99- 


HEATER  CYCLE:  PRESSURE  RATIO2  4.79 
TOTAL  PRESSURE  DROP2. 9 3 
AIR  FLOW ( L'om/ s )  =  39-22 
FUEL  FLOW ( Ltd/ s )  =  3*25 

STATE  TEM?( deg  R) 

lh  535. 0 

2h  835.2 

3h  129-2.5 

4h  4647-3 

5h  1321.1 

6h  945.0 

7h  760.0  ,  , 

COMPONENT  DP?  DKD  VCL(cu  - 
RSCUP  .235  *850  2.2 

TURBINE 
HEATER 
COMPRESSOR 


.098  .694 


VCL( cu  ft) 

2.2  SF?=0.84 

0.3  E?F( ?) =0 . 89 
5.0 

11.9  E??(?) =0.35 


RANKINS  CYCLE:  SATURATION  PRESSURE2  199.9 
STEAM  FLOW ( Ldm/ s ) =  10.76 
FLOW  RATI0*0 .06l3 


STATE 

Is 

2s 

3? 

3s 

4s 

COMPONENT 
BO  Ix«:R 
TURBINE 
r.riNnPNSZR 


TEM?( deg  R) 
540.0 
555-0 
341.7 
1103.1 
54o.o 
VOL(  cu  fJ 
1330.6 
34.0  1 

2<?7.1 


PRESSURE! os i) 

0.5 
199*9 
199.9 
199.9 
0.5  ac 


ZFF! IS) =0.38 


QUAL=C.375 


